Molecular Hydrogen Emission from Protoplanetary Disks II. 
Effects of X-ray Irradiation and Dust Evolution 
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ABSTRACT 

Detailed models for the density and temperature profiles of gas and dust in protoplanetary 
disks are constructed by taking into account X-ray and ultraviolet (UV) irradiation from a central 
T Tauri star, as well as dust size growth and settling toward the disk midplane. The spatial and 
size distributions of dust grains in the disks are numerically computed by solving the coagulation 
equation for settling dust particles, with the result that the mass and total surface area of dust 
grains per unit volume of the gas in protoplanetary disks are very small, except at the disk 
midplane. The level populations and line emission of molecular hydrogen are calculated using 
the derived physical structure of the disks. X-ray irradiation is the dominant heating source of 
the gas in the inner disk region and in the surface layer, while the far UV heating dominates 
otherwise. If the central star has strong X-ray and weak UV radiation, the level populations 
of molecular hydrogen are controlled by X-ray pumping, and the X-ray induced transition lines 
could be observable. If the UV irradiation is strong, the level populations are controlled by 
thermal collisions or UV pumping, depending on the properties of the dust grains in the disks. 
As the dust particles evolve in the disks, the gas temperature at the disk surface drops because the 
grain photoelectric heating becomes less efficient, while the UV radiation fields become stronger 
due to the decrease of grain opacity. This makes the level populations of molecular hydrogen 
change from local thermodynamic equilibrium (LTE) to non-LTE distributions, which results in 
changes to the line ratios of molecular hydrogen emission. Our results suggest that dust evolution 
in protoplanetary disks could be observable through the line ratios of molecular hydrogen. The 
emission lines are strong from disks irradiated by strong UV and X-rays and possessing small 
dust grains; such disks will be good targets in which to observe molecular hydrogen emission. 



Subject headings: line: formation 
protoplanetary disks 

Introduction 



molecular processes — radiative transfer — planetary systems: 



Observations of thermal dust continuum emis- 
sion (e.g., Kenyon & Hartmann 1987; Beckwith & 
Sargent 1993; Kitamura et al. 2002) and, more di- 
rectly, images of light scattered from dust grains 
(e.g., Roddier et al. 1996; Itoh et al. 2003; Duch- 
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ene et al. 2004) have revealed that young stellar 
objects have circumstellar disks. In addition, a 
gaseous component has been detected from the cir- 
cumstellar disks around T Tauri stars (e.g., Carr 
1989; Dutrey et al. 1997, 2007; Thi et al. 2004; 
Qi et al. 2006; Najita et al. 2007; Bergin et 
al. 2007). Recent high spectral resolution and 
high sensitivity observations have made it possi- 
ble to detect line emission of molecular hydrogen 
gas, which is the major component of the gas in 
the protoplanetary disks (Thi et al. 1999, 2001a, 
b; Richter et al. 2002; Sheret et al. 2003; Sako 
et al. 2005; Weintraub et al. 2000; Bary et al. 
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2002, 2003; Itoh et al. 2003; Herczeg et al. 2002, 
2004, 2006; Bergin et al. 2004). Furthermore, de- 
tection of molecular hydrogen lines in the near- 
and mid-infrared wavelength bands has been re- 
ported towards T Tauri stars and Herbig Ae/Be 
stars very recently (Weinstraub et al. 2005; Bary 
et al. 2007, in preparation; Richter et al. 2007; 
Bitner et al. 2007, in preparation). In a previous 
paper (Nomura & Millar 2005, hereafter Paper I), 
we constructed a model for molecular hydrogen 
emission from a protoplanetary disk that is irra- 
diated by strong ultraviolet (UV) radiation from 
a central star, and whose dust component has the 
same properties as dense molecular cloud dust. 

Now, it is known observationally that many 
young stellar objects emit strong X-ray radiation 
(Koyama et al. 1994; Feigelson & Montmerle 1999; 
Tsujimoto et al. 2002; Imanishi et al. 2003; Get- 
man et al. 2005; Preibisch et al. 2005) which 
ionizes hydrogen gas and could be one of the im- 
portant heating sources of gas in protoplanetary 
disks, in addition to the grain photoelectric heat- 
ing induced by far UV radiation from the central 
star (cf. Kamp & Dullemond 2004; Jonkheid et al. 
2004; Dullemond et al. 2007; Paper I). Actually, 
some model calculations show that X-ray irradia- 
tion can heat the gas to very high temperatures at 
the surface layer of the disks (Glassgold & Najita 
2001; Gorti & Hollenbach 2004; Glassgold et al. 
2004; Kamp et al. 2005). Furthermore, it has been 
suggested that fast secondary electrons produced 
by X-ray ionization, similar to the electrons in- 
duced by cosmic-ray ionization, can pump molec- 
ular hydrogen into excited electronic states and 
may be important, for example, in extra-galactic 
objects, Herbig-Haro objects, and supernova rem- 
nants (e.g., Shemansky et al. 1985; Gredel et al. 
1989; Gredel & Dalgarno 1995; Tine et al. 1997). 
This X-ray pumping could be observable toward 
disks which are irradiated by strong X-ray radi- 
ation from their central stars (e.g., Bergin et al. 
2004). 

As disks evolve, it is believed that the dust par- 
ticles in the disks coagulate and settle toward the 
disk midplane as the first step of planet formation. 
The dust dynamics in this stage have been studied 
theoretically in many works (e.g., Wcidcnschilling 
1980, 1997; Nakagawaet al. 1981, 1986; Mizuno et 
al. 1988; Mizuno 1989; Cuzzi et al. 1993; Schmitt 
et al. 1997; Nomura & Nakagawa 2006). This 



dust evolution is expected to affect observational 
properties of the disks, and some model calcula- 
tions have been done in order to study the effect 
(e.g., Miyake & Nakagawa 1993, 1995; DAlessio 
et al. 2001, 2006; Dullemond & Dominik 2004; 
Jonkheid et al. 2004, 2006, 2007; Retting et al. 
2006; Aikawa & Nomura 2006). In addition, some 
numerical calculations of the dust evolution have 
been done by solving the coagulation equation for 
settling dust particles, and its effect on the spec- 
tral energy distribution of thermal dust emission 
disks investigated (Suttner & Yorke 2001; Tanaka 
et al. 2005; Dullemond & Dominik 2005). In this 
paper, we further examine the effects of the dust 
evolution on the physical structure of the gas in 
the disks and on molecular hydrogen emission by 
using both a simple dust model and a numerical 
calculation of the coagulation equation. 

Historically, line emission from molecular hy- 
drogen has been observed towards various kinds of 
astronomical objects, such as shock fronts associ- 
ated with star forming regions, reflection nebulae, 
planetary nebulae, supernova remnants, external 
galaxies, and so on (e.g., Beckwith et al. 1978; 
Brown et al. 1983; Hasegawa et al. 1987; Bur- 
ton et al. 1992). The observed line ratios probe 
the physical properties of these objects as they 
reflect the excitation mechanisms of molecular hy- 
drogen, e.g., thermal collisions, ultraviolet and X- 
ray pumping, and formation pumping (e.g., Black 
& van Dishoeck 1987; Sternberg & Dalgalno 1989; 
Tanaka et al. 1989; Tine et al. 1997; Takahashi 
& Uehara 2001). In this work we propose a possi- 
ble observational diagnostic of the dust evolution 
in protoplanetary disks using line spectra and the 
line ratios of molecular hydrogen. 

In the following sections, we model the density 
and temperature profiles of the gas and dust in 
protoplanetary disks, taking into account the X- 
ray and UV irradiation from a central star, as well 
as dust growth and settling towards the disk mid- 
plane. Then, using the physical structure, we cal- 
culate the level populations and line emission of 
molecular hydrogen. In §2, we introduce the mod- 
els we use in this work. For the dust evolution, 
we use both a simple model and a more realistic 
model in which we solve the coagulation equation 
for settling dust particles ( £12. ip . The X-ray and 
UV radiation fields are computed 1+1 dimension- 
ally ( §2.2j) . the density and temperature profiles 
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of the gas and dust in the disks are obtained by 
assuming vertical hydrostatic equilibrium and lo- 
cal thermal and radiative equilibrium ( §2.3p . The 
level populations of molecular hydrogen are calcu- 
lated under an assumption of statistical equilib- 
rium from which we get the molecular hydrogen 
emission by solving the radiative transfer equa- 
tion ( §2.4|) . In §3, we present the resulting dust 
size and spatial distributions ( §3.1j) . the physical 
structure of the disks f H3.2p . the level populations 
of molecular hydrogen ( £|3.3[) . and the line spectra 
and line ratios of molecular hydrogen ( §3.4p . in 
which the effects of the X-ray irradiation and the 
dust evolution are discussed. Finally, the results 
are summarized in §4. 

2. Models 

2.1. Spatial and Size Distributions of Dust 
Particles 

The physical structure and the molecular hy- 
drogen emission of the disks are affected by the 
dust model in various ways; for example, the UV 
radiation field through dust extinction ( §2.2|) . 
the dust temperature through optical proper- 
ties of dust grains, the gas temperature through 
grain photoelectric heating and dust-gas collisions 
( ^2.31) . and the molecular hydrogen formation rate 
on dust grains ( §2.4p . In this paper we use the 
following two types of model for the dust size 
and spatial distributions. In model A we adopt 
a very simple assumption in order to understand 
the basic properties of the effects of dust evolu- 
tion. In model B we consider a more realistic case 
by numerically solving the coagulation equation 
for settling dust particles. 

In both models the shape of the dust particles is 
simply assumed to be a compact sphere (see e.g., 
Kozasa et al. 1992; Ossenkopf 1993; Ormel et al. 
2007 for fractal dust aggregate models). Making 
use of the resulting spatial and size distributions 
of dust grains, the dust absorption (k„) and scat- 
tering (av) coefficients at each position of the disk 
are computed by means of the Mie theory (Bohrcn 
& Huffman 1983), with the dust particles assumed 
to consist of silicate, carboneous grains, and water 
ice (see Paper I for details). 



2.1.1. Model A 

In model A we assume that the dust grains have 
a spatially uniform distribution and the mass frac- 
tional abundance of the dust with respect to the 
gas is fixed at each position in the disk (i.e. the 
dust grains are well mixed with the gas). The 
size distribution is set to be dn/da tx a~ 35 (a is 
the radius of dust grain) with the maximum grain 
radii of a max = 10/im, 1mm, and 10cm. The min- 
imum radius is set to be 0.01/xm for all models. 
The amount of small dust grains decreases with in- 
creasing maximum grain radii as we keep the mass 
fractional abundance of the dust grains to the gas 
fixed (see also §2.1.3) (e.g., Miyake & Nakagawa 
1993; DAlessio et al. 2001; Aikawa & Nomura 
2006). 

2.1.2. Model B 

In this model the spatial and size distributions 
of dust grains are obtained by solving coagulation 
equations for various sizes of settling dust parti- 
cles, 

^ + J^M*)] = -m^)X>(*, iMi) 

3=1 

1 

where ip(i) is the mass density of dust particles 
in a mass bin i whose summation is equal to the 
total mass density of the dust particles at a given 
position and time as 

n 

PdustQz, z,t) = y] <p{x, z, t, i). (2) 

i=l 

Here we briefly summarize the dust evolution 
model; more details can be found in Nomura & 
Nakagawa (2006). Now, V z (i) in equation JT]) is 
the vertical velocity and vtii is the typical mass of 
a particle in a mass bin i. The symbol (3(i,j) is 
related to the sticking rate of two colliding dust 
particles, given by 

/3(i,j) = 7r(<Zi + ajfSVps/mirrij, (3) 

where aj is the radius of a dust particle in a mass 
bin i, and we simply assume the sticking proba- 
bility of p s = 1 in this paper. The fragmentation 
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of dust particles is simply neglected in this work. 
For the relative velocity between two colliding par- 
ticles, SV, we adopt 

SV = (SVg + SV 2 + SV 2 + 5V% ) 1/2 . (4) 

The symbol SV B = {8kT d /Tr) 1 / 2 (l/m, l + 1/mj) 1 / 2 
{k is Boltzmann's constant and Td is the dust tem- 
perature) is the relative velocity caused by the 
thermal Brownian motion. The symbols SV Z = 
V z (i) ~ V z (j) and SV X = V x (i) - V x (j) (V x (i) is the 
local radial velocity component of dust particles 
with mass TOj, arising from angular momentum 
loss via gas-dust friction) are the velocity differ- 
ences in the vertical and radial directions, respec- 
tively. Finally, SVr is the turbulence-induced rel- 
ative velocity (see Nomura & Nakagawa 2006 for 
more details). The disk is simply assumed to be 
completely quiescent or turbulent, with SVr = 
in a quiescent disk model. 

The mass flux in equation ([T]) is given by 

K,(i) ¥ ,(i) = -^(i) (5) 
Ap 

in a quiescent disk, where Qk is the Keplcrian fre- 
quency and p is the gas density. For the drag co- 
efficient between the gas and dust particles, A, we 
adopt A — c s /p s a for a < l g and A — 3c s l g /2p s a 2 
for a > l g , following Epstein's and Stokes' laws, 
respectively, where c s , p s , and l g are the sound 
speed of the gas, the solid density of a dust par- 
ticle, and the mean free path of gas particles, re- 
spectively. The mean velocity of the dust parti- 
cles in the vertical direction, V z (i) = tt^z/Ap, is 
obtained by balancing the gas-dust friction force, 
ApV z (i), and the gravitational force in the vertical 
direction, tt^z. In a turbulent disk, the mass flux 
is written as 

V z (i)v(i) = -%(i)-D p d J£m, (6) 
Ap oz 

following the gradient diffusion hypothesis. For 
the turbulent diffusivity, we adopt Dq — a'c s H/(l + 
flx/Ap), where H — c s0 /£Ik (c s0 is the sound 
speed at the disk midplane) is the scale height of 
the disk and we set a' — 10~ 4 in this paper. 

The global radial motion of the dust particles 
toward the central star is not taken into account in 
this work. At the disk surface this simplified treat- 
ment is applicable because large dust particles set- 
tle toward the disk midplane more rapidly than 



they move toward the central star (Figs. [3] and [4] 
in M3 . 1 1 actually show that large particles cannot 
stay in the surface layer). We note that the radial 
motion is negligible for the small dust particles 
which couple with the gas efficiently via friction 
force (e.g., Adachi et al. 1976; Weidenschilling 
1977; Takeuchi & Lin 2005). Therefore, neglect of 
radial global motion will not affect our results of 
molecular hydrogen lines which are mainly emitted 
in the disk surface (see Paper I). In a completely 
quiescent disk the radial motion is negligible all 
over the disk (e.g., Nakagawa et al. 1986). Now, 
in order to avoid making very large particles which 
do not couple with the gas and should fall on to- 
wards the central star, we remove particles larger 
than some critical radius by simply assuming that 
when the dust particles grow large enough so that 
they cannot be trapped in a turbulent eddy, they 
gain very rapid radial motion. The critical radius 
is estimated as a cr ;t = c s p gas / p s Qk for a < l g and 
a C rit = (3c s p gas ^/2p s f2 K ) 1/2 for a > l g by com- 
paring the friction time between the gas and dust 
particles, r/ = l/Ap, and the turnover time of the 
largest turbulent eddy, r c ddy = f/^K- We note 
that the dust particles grow to be larger than a cr it 
only close to the midplane of a turbulent disk (see 

As the initial condition of the calculations, we 
set the dust particles to be well-mixed with the gas 
and the dust-to-gas mass ratio to be spatially uni- 
form. For the initial size distribution of the dust 
grains, we adopt dust model A with a max = 10/im, 
which is similar to the dust model of dense molec- 
ular clouds (e.g., Weingartner & Draine 2001). 
We also assume that the disk is surrounded by 
a dense molecular cloud with gas number density 
"-out = 10 4 cm -3 , and consider a continuous input 
of dust particles from the molecular clouds to the 
disk as a boundary condition at the disk surface 
(z = z coa g, see below and §3.ip for the calculations 
of dust evolution (Eq. [T]). Dust particles fall on 
to the disk because the pressure gradient force is 
negligibly small for them (e.g., Landau et al. 1967) 
and cannot sustain them against the gravitational 
force of the central star (we note that the gas is 
assumed to be stationary due to the pressure gra- 
dient force) . This input of dust particles has great 
influence on structure of the disk surface and ob- 
servational properties of the disk (see [J3|) . The to- 
tal mass of dust particles infalling from the cloud 
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to the disk in 10 6 yr (the time for which the calcu- 
lations are performed) is ~ 5 x 10 _5 M© , which cor- 
responds to a dust mass in a spherical cloud with a 
radius of 2,500 AU and a gas density of 10 4 cm~ 3 , 
and 1/3 of the initial mass of dust grains in the 
disk. In the calculations for the dust coagulation 
and settling we set 23 radial grids logarithmically 
for x = 0.2 — 100 AU, and 50 vertical grids for the 
region where the coagulation becomes significant 
[z < z coag , see i )3.ip . The dust coagulation equa- 
tion is solved using fixed density and temperature 
profiles of dust and gas because self-consistent cal- 
culations coupled to the time evolution of the dust 
and gas profiles are very time-consuming. Instead, 
we get the dust distribution and the temperature 
and density profiles by iterating the calculations 
only once under an assumption that the tempera- 
ture and density profiles turn into an equilibrium 
state very quickly; that is, we (a) calculate the 
dust gas profiles using the initial dust distribution, 

(b) solve the coagulation equation using the tem- 
perature and density profile in (a), (c) compute 
the dust and gas profiles again using the dust dis- 
tribution in (b), (d) obtain the dust distribution 
by solving the coagulation equation with the tem- 
perature and density profiles in (c) , and (e) finally 
get the dust and gas profiles using the dust distri- 
bution in (d). In the calculations in [J3]we com- 
pare the gas temperature and density profiles in 

(c) and (e), and checked that the errors are within 
30% for the gas temperatures and 80% at most 
for the gas densities in both quiescent and turbu- 
lent disks. These errors are relatively large at the 
upper layer of the disks where the molecular hy- 
drogen is photodissociated via UV radiation from 
the central star, so the errors in molecular hydro- 
gen emission from the disks (see ^3.4p are smaller 
and within 15% for all lines. 

2.1.3. Parameter for Total Surface Area of Dust 
Particles, /dust 

We introduce a parameter, /dust, that repre- 
sents the total surface area of dust particles per 
unit volume of the gas at each position in the disk 

(x,z), 

/dust (a, z) = A tot (x, z)/A tot ,o, (7) 

where 

A tot (x,z) = U-Ka 2 ^§^-da. (8) 



The size distributions dn/da oc ip(i)/af in dust 
model B, and A to t,o is calculated using the dust 
model A with a max = 10/im, which is similar to 
the dust model of dense molecular clouds and is 
used as the initial condition for dust model B. 
This parameter controls the physical disk struc- 
ture and molecular hydrogen emission from the 
disks because the grain opacity, the grain pho- 
toelectric heating rate, and the formation rate of 
molecular hydrogen on grain surface are roughly 
proportional to it. Thus, the dust and gas temper- 
atures, the UV radiation field, and the abundance 
of molecular hydrogen are related to /dust (see the 
following subsections). In dust model A, the pa- 
rameter /dust has values of 1, 0.1, and 0.01 for the 
maximum dust size of a max = 10/im, 1mm, and 
10cm, respectively. In model B, /d us t decreases 
with dust size growth and settling (except at the 
disk midplane), which leads to a decrease of the 
grain photoelectric heating rate, and thus the gas 
temperature and so on (see fc|3.1l for more details). 

2.2. X-ray and Ultraviolet Radiation Fields 

Observations have shown that many T Tauri 
stars emit strong X-ray (e.g., Koyama et al. 1994; 
Feigelson & Montmerle 1999; Tsujimoto et al. 
2002; Imanishi et al. 2003; Getman et al. 2005; 
Preibisch et al. 2005) as well as strong ultravio- 
let (UV) radiation (e.g., Herbig & Goodrich 1986; 
Herbst et al. 1994; Valenti et al. 2000). For 
the X-ray radiation from the central star we use 
a model which reproduces observational data to- 
ward a classical T Tauri star, TW Hydrae (cf. 
Kastner et al. 2002; Stelzer & Schmitt 2004). 
Retrieving the archived XMM-Newton data, we 
fit the spectrum with a two-temperature thin- 
thermal plasma model (mckal model; Mewe et 
al. 1985; Kaastra et al. 1992; Licdahl et al. 
1995) which is often used in order to reproduce 
observed X-ray spectrum of T Tauri stars. The 
derived best-fit parameters are kT% = 0.8keV and 
kT-i — 0.2keV for the plasma temperatures, and 
Vh = 2.7 x 10 20 cm~ 2 for the foreground inter- 
stellar hydrogen column density. The total X- 
ray luminosity of the spectrum corresponds to 
Lx ~ 10 30 erg s _1 . In Figure [1] the resulting 
model spectra is plotted. The adopted stellar UV 
radiation field model is also based on observations 
towards TW Hydrae and analyses by Herbst et al. 
(1994), Costa et al. (2000), Bergin et al. (2003), 
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Fig. 1 . — The model spectra of the X-ray radiation 
at the central star, which reproduces observation 
toward a classical T Tauri star, TW Hya (d = 
56pc). 

Herczeg et al. (2002), and Ardila et al. (2002). 
The model consists of photospheric black body ra- 
diation, hydrogenic thermal bremsstrahlung radi- 
ation, and strong Ly a line emission (see Appendix 
C of Paper I). The total FUV (6eV < hv < 13eV) 
luminosity corresponds to Lfuv ~ 10 31 erg s _1 . 
The interstellar UV radiation field is taken into 
account, but its contribution is negligible under 
the strong UV irradiation from the central star 
(see Paper I for details of the UV radiation field 
in the disk). We note that although we use the 
X-ray and UV radiation of TW Hya (as it is one 
of the most well-observed T Tauri stars) , our disk 
model is to be more widely applicable. 

The X-ray and UV fields in the disk are calcu- 
lated in 1+1 dimensions in the radial and vertical 
directions (see also Paper I) as 

f R 

F U! r(R, 6) = /F 1 , !8tar exp(-r I/ ,fl), t u ,r= / XvpdR, 

Jr. 



(9) 



and 



F„ tg (x, z) = F^isrf exp[-r 1 , iZ (z 00 )] 
+2tt / °° a v {x,z')p{x,z')F v , R (x,z')e- T ^ z 'Uz' , 



Tv.z(z') = / X v pdz" ', 



(10) 



where the direct radiation from the central star 
is calculated fully, while a plane-parallel approx- 
imation in the vertical direction is adopted for 



the calculation of the scattering process, which 
could result in overestimating the radiation fields 
in the disks. In the above equations -F^star is the 
specific radiation field at the stellar surface and 
/ = (R*/R) 2 accounts for the geometrical dilution 
of the radiation field. F^isrf is the FUV interstel- 
lar radiation field, and we set Fx r ay, isrf = in 
the calculation of the X-ray radiation field. t u ,r 
and Ty jZ are the specific optical depths from the 
stellar surface (i?„, 0) to a point (R, 6) and from 
a point (x, z) to (x,z'), respectively, p is the gas 
density, and Xv is the monochromatic extinction 
coefficient defined by the absorption (k„) and scat- 
tering (ov) coefficients as Xv = K v + Vv In order 
to treat X-ray extinction, we adopt the attenu- 
ation cross section at an energy E of <T at t(E) = 
Cph(E)+<jCom(E), where cr p h is the total photoion- 
ization cross section due to all elements per hydro- 
gen nucleus and acom is the incoherent Compton 
scattering cross section of hydrogen. For the cross 
section <j p h we adopt a broken power-law model 
given by Maloney et al. (1996; see also Wilms et 
al. 2000), and ocom is calculated based on Mc- 
Master et al. (1969, http://cars9.uchicago.edu/~ 
newville/mcbook/). In calculating equations (|9|) 
and (TTU]), we adopt xx ray = e&tt (E)/m p and 
cxray = <Jcom (F) / m p , where m p is the proton 
mass. 



2.3. Physical Structure of the Disks 

We model an axisymmetric disk surrounding a 
central star with the physical parameters of typical 
T Tauri stars; a mass of M* = 0.5Af Q , a radius of 
i?» = 2_R©, and a temperature of T* = 4000K 
(e.g., Kenyon & Hartmann 1995). 

The gas temperature and density distributions 
of the disk are obtained self-consistently by iter- 
atively solving the equations for hydrostatic equi- 
librium in the vertical direction and local thermal 
balance between heating and cooling of gas (see 
Paper I for details). The vertical hydrostatic equi- 
librium is represented by an equation, 



dP 

dz 



GAL 



-P 



(x 2 + z 2 ) 3/2 ' 



(11) 



G is the gravitational constant, and P is the 
gas pressure given by P — pkT/m^, where p, 
T, k, and m M are the density and tempera- 
ture of the gas, Boltzmann's constant, and the 



6 



mean molecular mass, respectively. The condi- 
tion, Jj^ p(x,z)dz = X(x), is imposed, where 
we set ,0(2;, = 5.0 x 1CT 19 g cm" 3 (n tot ~ 
3 x 10 5 cm~ 3 ) as the boundary condition. The 
surface density at a disk radius x, S(x), is defined 
by assuming a constantly accreting viscous disk 
model and equating the gravitational energy re- 
lease of accreting mass to the thermal heating via 
viscous dissipation at the disk midplane, 



-Sac S gO K = 



3GM*M 

47TX 3 



X 



1/2 



(12) 



where c s0 and f2 K = (GM^/x 3 ) 1 / 2 represent the 
sound speed at the midplane and the Keplerian 
frequency, respectively. A viscous parameter of 
a = 0.01 and a constant mass accretion rate of 
M = 10 -8 Mq yr _1 are adopted here. 

The gas temperature, T, is obtained by assum- 
ing detailed energy balance at each position in the 
disk, 



FUV 



+ r 



Xray 



L„ T + A liE 



0. 



(13) 



where we include grain photoelectric heating in- 
duced by FUV photons, Tpuv, X-ray heating 
caused by hydrogen ionization, Txray, gas-grain 
collisions, A gr , and radiative cooling by line tran- 
sitions, Aii ne for the gas heating and cooling pro- 
cesses. The X-ray heating rate, Txray, is calcu- 
lated as 

Txray = n tot f h Hx, (14) 

where n to t is the total number density of hydrogen 
nuclei, and fh is the heating efficiency, namely, the 
fraction of absorbed energy that goes into heating 
the gas. We adopt fh = 0.1 for atomic hydrogen 
and fh — 0.4 for molecular hydrogen (Maloney 
et al. 1996; Gorti & Hollenbach 2004). H x is 
the local X-ray energy deposition rate per particle, 
given by 



H x = 



a ph (E)F x (E)dE, 



(15) 



where a p h(E) is the total photoionization cross 
section due to all elements per hydrogen nucleus 
at energy E. The symbol Fx(E) is the X-ray 
energy flux at each position in the disk, and 
E m i n = O.lkeV and E maK — lOkeV are adopted 
for the minimum and maximum energy (see Fig.[T] 
in £ |2.2p . We note that the viscous heating is not 



taken into account in the energy balance (Eq. [13] ) 
because it is not dominant (at the disk surface) if 
a = 0.01 (Glassgold et al. 2004). 

For radiative cooling by line transitions, we con- 
sider the Ly a transition of atomic hydrogen and 
the metastable transition of OI (A 6300A) in addi- 
tion to the fine-structure transitions of OI (63/im) 
and CII (158^m), and the rotational transitions 
of CO. In order to calculate the Ly a line cool- 
ing, we make use of the table of level populations 
of atomic hydrogen for various electron densities 
and temperatures, given by Storey & Hummer 
(1995). The collisional de-excitation rate coeffi- 
cient is taken from Hollenbach & McKee (1989) 
for calculation of OI A 6300 line cooling. Paper 
I gives details of calculations of the OI and CII 
fine-structure, and CO rotational transition line 
cooling. 

The spatial and size distributions of dust grains 
affect the gas temperature through the grain pho- 
toelectric heating, Tpuv, and the energy exchange 
between gas and dust particles through collisions, 
A gr . Both rates are roughly proportional to the 
parameter which represents the total surface area 
of the dust particles, /dust, given in t|2.1.31 In 
this paper we simply set r FUV = /dustFpuv.o 
and A gr = /dustAgr,o, where the heating/cooling 
rates with subscript '0' are calculated by using 
the models given in Paper I in which we used the 
dense cloud dust model (see also Aikawa & No- 
mura 2006). 

The dust temperature profile is important for 
determining the disk structure because the gas 
temperature is well coupled to the dust temper- 
ature in the dense region near the midplane of the 
disks. We obtain the dust temperature by assum- 
ing local radiative equilibrium between absorption 
and reemission of radiation by dust grains at each 
position in the disk. The intensity is calculated by 
solving the axisymmctric two-dimensional radia- 
tive transfer equation by means of the short char- 
acteristic method in spherical coordinates (Dulle- 
mond & Turolla 2000; Nomura 2002). As heating 
sources, we consider the radiative flux produced 
by the viscous dissipation (a- viscous model) at the 
disk midplane, and the irradiation from the central 
star (see Paper I for details). The dust evolution 
in the disks affects the dust temperature through 
the change in grain opacity ( §2.1|) . 
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2.4. Level Populations and Line Emission 
of Molecular Hydrogen 

In order to obtain the molecular hydrogen emis- 
sion from the disk, we first calculate the abun- 
dance and the level populations of the elec- 
tronic state of molecular hydrogen in a statis- 
tical equilibrium state, based on Wagenblast & 
Hartquist (1988), as 



is not too strong, the gas temperature is moder- 
ately high, and the surface area is large; see Paper 
I), but a full calculation of X-ray pumping and the 
subsequent radiative cascade should be done in fu- 
ture. The symbol £x in equation (fTT|) is the total 
hydrogen ionization rate given by 



a ph {E)F x (E)dE, 



(18) 



ni(R 2 



E 



■ft. 



■llr, 



+fc o +H 2 n(0)n ; (H 2 ) 



(H 2 ) ( A m l+p m t+J m l+^2 nsC ™l ) +«(H)i?form 
^ s ' 

(16) 



where Ai m is the Einstein A-coefficient for sponta- 
neous emission from level I to level m and Cf m is 
the collisional transition rate with collision partner 
s. film represents the effective rate for transition 
I — > m via ultraviolet pumping followed by radia- 
tive cascade, and i?diss,z is the photodissociation 
rate of hydrogen molecules in level I. Rf orm ,i is the 
effective formation rate of H2 in level / on grain 
surfaces. The endothermic reaction O + H2 — > 
OH + H, which destroys molecular hydrogen in 
high temperature regions, is also taken into ac- 
count (see Paper I for details). In addition, we 
consider the effective transition rate, 7; m , via X- 
ray pumping of molecular hydrogen followed by 
radiative cascade. X-ray irradiation from the cen- 
tral star ionizes the gas to produce photoelectrons 
and subsequently secondary electrons in the disk. 
Through collisions they excite molecular hydrogen 
to singlet and triplet electronic states, followed by 
radiative cascade down into the ground electronic 
state (e.g., Gredel & Dalgarno 1995; Tine et al. 
1997; Bergin et al. 2004). In this paper we simply 
use the entry efficiency otj^v, J) from the levels 
(v = 0, Ji) to (v, J) for fractional ionization of 
10~ 4 , tabulated in Tine et al. (1997), in order to 
estimate the rate 7; TO as 



7fm = Cxaj,(Wm, Jm)- 



(17) 



This simplified treatment will not cause significant 
error at the middle layer of the outer disk where 
molecular hydrogen lines are mainly emitted (the 
line fluxes are strong where the UV radiation field 



where cr p h is the total photoionization cross section 
due to all elements per hydrogen nucleus, and N sec 
is the number of secondary ionizations of hydro- 
gen per unit energy produced by primary photo- 
electrons and we put N sec = 26/keV in this pa- 
per (e.g., Verner & Yakovlev 1995; Maloney et al. 
,*1996; Gorti & Hollenbach 2004). We note that the 
effect of interstellar cosmic-ray ionization is not 
taken into account in the calculation of the pump- 
ing process, but it will not affect the results as the 
X-ray ionization rate is much higher at the disk 
surface (see £13.2. 1|) . Possible reactions induced by 
X-rays are ignored, and a simple chemical network 
given in Wagenblast & Hartquist (1988) (plus the 
reactions O + H 2 -> OH + H and OH + hv -> O 
+ H) is adopted as in Paper I. This neglect will not 
affect the resulting molecular hydrogen abundance 
since the photodissociation by UV radiation from 
the central star or the above-mentioned reaction 
with atomic oxygen is more efficient for destroy- 
ing molecular hydrogen than the X-ray-induced 
photoionization and other related reactions in our 
model (see e.g., Maloney et al. 1996 for XDR 
chemistry, and also e.g., Aikawa & Herbst 1999, 
2001; Markwick et al. 2002 for more detailed disk 
chemistry including the X-ray photoprocess). 

The spatial and size distributions of dust grains 
affect the formation rate of molecular hydrogen. 
The rate is roughly proportional to the total sur- 
face area of the dust particles, /dust, given in 
12.1.31 so we simply set Rf olm j = /dust-Rform.i.o- 
Here we use the model in Paper I in order to calcu- 
late -Rf O rm,2,0- The total formation rate is given by 

Ei flfamj = 7-5 x 10- 18 f dust T°- 5 e H2 (T d )n tot n(H) 
cm~ 3 s^ 1 , where T is the gas temperature and 
£H 2 (2rf) is the recombination efficiency of atomic 
hydrogen on dust grains as a function of the dust 
temperature, T d (Cazaux k Tielens 2002, 2004; 
see also Pirronello et al. 1999; Zecho et al. 2002). 

Making use of the physical properties obtained 
in the previous subsections and the level popula- 
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tions, we calculate emission (from levels u to V) 
of molecular hydrogen from the disks by integrat- 
ing the radiative transfer equation (see Paper I for 
details), 

1 [ x °™ [ z °° 
F uI = ~a — 71 / ^xdx / fj u i{x, z)dz, (19) 

where fj u i(x,z) is the emissivity of the transition 
line at (x,z) times the effect of absorption in the 
upper disk layer, given by 

Vui(x,z) = n u (x, z)A u J^- exp(-r„ ; (x, z)). 

(20) 

t u i(x,z) is the optical depth from z to the disk 
surface z^, at the frequency v u i, 

T ui (x,z)=J^ Xui(x,z')dz', (21) 

where Xui is the total extinction coefficient, 

Xui = PXv ul + {niBiu - n u B u i)^ u i^^-. (22) 

In these equations, A u i and B u i are the Einstein 
coefficients, n u and n\ are the number densities of 
the upper and lower levels, respectively, and is 
the line profile function. The energy difference be- 
tween the levels u and I corresponds to hv u i. The 
symbol Xv n i IS ^ ne extinction coefficient of dust 
grains (see §2.11 and §2.2|) at the frequency v u i , 
and p is the gas density. Here, the disk is assumed 
to be face on an observer, and we use the distance 
to an object of d = 56 pc for calculating the inten- 
sity in order to compare it with the observations 
towards TW Hya. Extinction by a foreground in- 
terstellar dust grains is not taken into account in 
the calculations. 

3. Results 

3.1. Spatial and Size Distributions of Dust 
Particles 

For dust model B we obtain the spatial and size 
distributions of dust particles by solving the coag- 
ulation equations for various sizes of settling dust 
particles in a quiescent or turbulent disk f ^2.1[) . 
In Figure [2] we plot the resulting profiles of (a) 
the parameter representing the total surface area 
of dust grains, /dust (defined in Eq. [7] in §2.1.3|) . 
and (b) the total dust density pdust normalized 



by the initial value pdust,o (defined in Eq. [2] in 
£ 12 . 1 . 2|) in the vertical direction at the disk radii 
of x — 1,10, and 100AU. The initial dust den- 
sity is simply proportional to the gas density (the 
dust particles are well-mixed with the gas) and 
corresponds to roughly 1% of the gas mass den- 
sity, Pdust,o ~ O.Olp, in this model. The solid 
and dashed lines are the profiles in quiescent and 
turbulent disks, respectively. The calculations are 
performed for 10 6 yrs, comparable to the typical 
age of classical T Tauri stars. We note that at 
t ~ 10 6 yr the dust coagulation process and set- 
tling motion (input from the upstream and out- 
put to the downstream) are almost in equilibrium 
state at each position in the disk, and the spa- 
tial and size distributions of dust particles do not 
change with time except in the region very close to 
the disk midplane. Therefore, the spatial and size 
distributions of dust particles in the surface layer 
presented in this subsection are applicable to older 
star-disk systems as well. Figure [5] shows that the 
mass and total surface area of dust grains per unit 
volume of the gas are much smaller than the ini- 
tial values. In the disk surface (z > z coag ; see be- 
low) where the density of particles is low enough 
so that the dust particles settle before they grow, 
Pdust/pdust,o and /dust are small due to the set- 
tling of dust particles toward the disk midplane. 
In the upper surface of the disk (z > Zf r i C ; see be- 
low), where the density is low enough that the 
gas friction force does not affect the motion of 
dust particles, the particles settle in the vertical 
direction with the free-fall velocity, V z = Vgf = 
[2GM*/(x 2 + z 2 ) 1 / 2 } 1 / 2 . In this region the nor- 
malized dust density /Odust//Odust,o (and the param- 
eter /dust) drop with decreasing z, inversely pro- 
portional to the gas (or the initial dust) density. 
Here, the dust particles are assumed to be con- 
tinuously falling on from the surrounding molec- 
ular cloud to the disk due to the gravitational 
force of the central star with a constant (time- 
independent) mass flux (= n ou tVff) (see fl2.1.2[l . 
At smaller z (z < zmc) where the gas density be- 
comes higher and the gas friction force controls the 
dust motion, the vertical velocity of dust particles 
becomes V z = il^z/Ap ( §2.1. 2| . and the normal- 
ized dust density /Cdust/pdust.o (and the parameter 
/dust) do not change very much in this region. The 
velocity changes from the free-fall velocity (Vs) 
to the terminal velocity (V z — ti^z/Ap) around 
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z = Zf ric = 0.5, 7.5, and 75 AU at the disk radii 
of a; = 1, 10, and 100AU, respectively, in this 
model. At even smaller z (z < z coag ) where the 
density is much higher and the collisional cross sec- 
tion becomes high enough for the dust particles to 
grow, the parameter /dust drops with decreasing 
z (and increasing density) because small particles 
disappear as a result of coagulation, while the nor- 
malized dust density Pdust/pdust.o does not change 
very much and increases close to the disk mid- 
plane due to the settling of the particles. Most 
of the dust mass settles at the disk midplane and 
Pdust/pdust,o > 1 at z ~ (not shown in this fig- 
ure). Throughout the calculations the total dust 
mass in the disk is equal to the initial dust mass 
plus the mass infalled from the cloud (minus the 
mass of particles with a > a cr j t removed near 
the midplane of the turbulent disk; see below). 
The difference between the quiescent and turbu- 
lent disks shows up most clearly in the parameter 
/dust at small z (z < z coag ; where the dust co- 
agulation is efficient) because the collision rate is 
higher in the turbulent disk owing to the turbu- 
lent induced relative velocity between the parti- 
cles, 5Vt (see §2.1.2p . The coagulation becomes 
efficient around z = z CO ag ~ 0.15 (0.2), 3.5 (4.0), 
and 65 (65) AU for the quiescent (turbulent) disk 
at x — 1, 10, and 100AU, respectively, in this 
model. 

In Figures [3] and H] we plot the resulting size dis- 
tributions of mass density of dust particles, <p(i), 
normalized by the initial dust density /Jdust,o, in 
quiescent and turbulent disks, respectively. Each 
figure shows the size distributions at (a) x = 1AU, 
t = 1 x 10 2 yr, (b) x = 1AU, t = 1 x 10 6 yr; 
(c) x = 10AU, t = 3 x 10 3 yr, (d) x = 10AU, 
t = 1 x 10 6 yr; (e) x = 100AU, t = 3 x 10 4 yr, 
and (f) x = 100AU, t = 1 x 10 6 yr. The time used 
in Figure a, c, and e is around the time when the 
size of the largest dust particles at the disk height 
of z ~ H becomes maximum in the quiescent disk 
model. The dot-dashed, dashed, and solid lines 
in each figure represent the size distributions at 
z ~ 2 coag , z ~ 2H, and z ~ H, respectively. In 
Figure SJd, d, and f, we also plot the size distribu- 
tions at z ~ 0.25-ff in thin solid lines. The disk 
scale heights are H = 0.044 (0.047), 0.51 (0.60), 
and 11 (11) AU for the quiescent (turbulent) disk 
at x = 1, 10, and 100AU, respectively. The thin 
dotted lines show the distribution of the initial 
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Fig. 2. — The vertical profiles of (a) the parame- 
ter for the total surface area of dust grains, /dust, 
and (b) the total dust density, pdust, normalized 
by the initial value in quiescent (solid lines) and 
turbulent (dashed lines) disks at the disk radii of 
x = 1, 10 and 100AU at 10 6 yr after the calcula- 
tions start. At the disk surface Pdust//?dust,o and 
/dust are small due to the dust settling towards 
the disk midplane. Near the disk midplane /dust 
is further smaller as a result of the dust coagula- 
tion, while Pdust/pdust.o increases due to the dust 
settling. 



10 



10" 

io- 2 
° io 4 

w 

a. 

S 10~ 6 

icr 8 

10 10 

10" 
10 2 

° 10- 4 



§ 10~ 6 



1 1 ,l I I I, I . ll I II, 
. (a)x=1AU,t=1x10 2 yr 


z=0.044AU 

Z-0.088AU 

z^D.15AU 

\ initial 






. ! . .1 . . .i . . .i !. . .i 


. i . .i . . Ii . . .i . . . 



dust radius [urn] 



10"' 
10 1t 

10° 

1(P 



(c) x=1 OAU, t=3x1 3 yr ^foAU 




dust radius [urn] 



° 10 4 >f 



^ 10 c 



i i-n — i i-n — i r-q— i r-q— i r-q— ■ 1-^—1 r-rf- 

(e) x=1 OOAU, t=3x1 4 yr ^=^1 AU — 

z=65AU - — 
initial 




10 1 10 2 10 3 
dust radius [urn] 



10" 
10 2 
°- 10 4 

Q. 

S 10~ 6 
10~ 8 
IO 10 

10" 
IO 2 

° 10- 4 



(b)x=1AU,t=1x10 6 yr 


z=0.044AU 

z=0.088AU 

z=0.15AU 

initial 












. .1 . . .1 . . .1 . . . 



dust radius [|im] 



10~ 6 



10"' 
10 1t 

10° 
IO 2 



(d)x=10AU,t=1x10 6 yr ^frjAU 

z=3^5AU 
initial 




dust radius [urn] 



1(f 



^ 10 c 
10* 
10" 1( 



(f)x=100AU,t=1x10V 




z=11AU 
z=22AU 
z=65AU 
initial 



"° 10 1 10 2 10 3 
dust radius [urn] 



10° 



Fig. 3. — The size distributions of mass density of dust particles, ip(i), normalized by the initial dust density 
Pdust,o at each disk radii, x, and time, t, in a quiescent disk. The dot-dashed, dashed, and solid lines 
represent the distributions at z ~ z coag , z ~ 2H, and z ~ H, respectively. The thin dotted lines show the 
initial distribution. At the disk surface the size distributions arc similar to those in molecular clouds, but 
the number density of the dust particles is much smaller than the initial value due to the dust settling. Near 
the disk midplane small particles disappear due to the dust coagulation and larger particles settle towards 
the midplane as time increases. 
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Fig. 4. — The same as Figure[3]but in a turbulent disk. The size distributions of dust particles at z ~ 0.25H 
are also plotted in thin solid lines in Figure b, d, and f. Near the disk midplane a certain number of large 
dust particles remain due to turbulent mixing. 
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condition (the dust model A with a max = 10/im). 
The figures show that at the surface layer above 
z ~ Zcoag the size distributions at t = 10 6 yr are 
similar to those in dense molecular clouds (the 
dust model A with a max = 10/im in this work) as 
the dust particles can not grow due to the small 
collisional rate, but the number density of the par- 
ticles is much smaller than the initial value due to 
the dust settling as mentioned above. We note 
that bumps of mass density of small dust parti- 
cles (a < 1/im) at z ~ z coag in early phases are 
remnants of the initial distribution. At smaller z 
(z < Zcoag), small dust particles disappear as they 
stick together to make larger particles. In the qui- 
escent disk the larger particles settle toward the 
disk midplane and disappear from the disk sur- 
face, z > H, as time goes on (Fig. [3]). Mean- 
while, in the turbulent disk a certain amount of 
large particles remain even at z > H at t — 10 6 yr 
(though most of them settle toward the midplane) 
because of the turbulent mixing which works so as 
to unify the size distributions in the vertical di- 
rection ( §2.1.21 Fig. H}. The cutoffs around the 
dust radii of a ~ 7 and 0.8 cm at the disk height 
of z ~ 0.25H in Figure [Uff and e correspond to 
the critical radii, a cr ;t, beyond which the particles 
cannot be trapped in a turbulent eddy and move 
toward the central star rapidly. The particles with 
a > a cr it are simply removed from the calculations 
(see glXU). 

3.2. Physical Properties of the Disks 

In this subsection we obtain the gas density 
and temperature distributions of the disk self- 
consistently by iteratively solving the equations 
for vertical hydrostatic equilibrium and local ther- 
mal balance between heating and cooling of gas 
( §2.3[) . The effects of the X-ray irradiation from 
the central star and the dust evolution on the 
physical properties of the disks are discussed in 
the following. 

3.2.1. Effect of X-rays 

First, in Figure [5] we plot the gas tempera- 
ture profiles in the vertical direction at the disk 
radii of x — 1, 10, and 100 AU, where the disk 
is irradiated by both of X-ray and UV radiation 
from the central star {solid lines). We also plot 
the profiles for a disk which is irradiated by X- 
ray radiation only {dashed lines) or UV radia- 



X-ray+UV 




Fig. 5. — The vertical temperature profiles of dust 
(thin dotted lines) and gas at the disk radii of x = 
1, 10 and 100AU for the irradiation models of X- 
rays + UV (solid lines), X-rays only (dashed lines), 
and UV only (dot-dashed lines). The dust model 
A with et max = 10/zm is used. The X-ray heating is 
dominant at the inner region and the very surface 
layer of the disk, while the FUV heating dominates 
in the middle layer and the outer region of the disk. 



tion only (dot-dashed lines) for comparison. The 
thin dotted lines are the dust temperature pro- 
files which are not affected by the UV or X-ray 
irradiation model. The dust model A with the 
maximum dust radius of a max = 10/im (see <j2.1[) 
is used throughout this sub-subsection. We note 
that the calculations are performed in the region 
where p > p(x, z^) = 5.0 x 10~ 19 g cm" 3 , and the 
position of z^ depends on the models (see §2.31) . 
The figure shows that the gas temperature is much 
higher than the dust temperature in the surface 
layer of the disk due to the X-ray and FUV heat- 
ing. The X-ray heating dominates the FUV heat- 
ing in the inner region and in the surface layer of 
the disk where direct irradiation from the central 
star is strong. Meanwhile, the FUV heating dom- 
inates the X-ray heating in the middle layer and 
in the outer disk. This is because the FUV radi- 
ation is scattered efficiently by dust grains, while 
the Compton scattering of X-ray radiation is inef- 
ficient in the energy range of E < IkeV (e.g., Igea 
& Glassgold 1999) in which T Tauri stars mainly 
emit X-rays (see Fig. [T|) . The gas temperature is 
almost the same as the dust temperature near the 
disk midplane where the density is high enough 
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that the gas and dust particles are well coupled 
through collisions. 

In Figure [6] we plot the vertical profiles of the 
heating and cooling rates at disk radii of (a) 1AU, 
(b) 10AU, and (c) 100AU, for a disk irradiated 
by both X-ray and UV radiation from the cen- 
tral star. The figures clearly show that the X-ray 
heating dominates in the inner region and in the 
surface layer of the disk, while the FUV heating 
dominates in the middle layer and the outer region 
of the disk. With regard to the cooling processes, 
radiative cooling dominates in the surface layer, 
while dust-gas collisions dominates near the mid- 
plane where the density is high. The main coolant 
at the surface layer changes as Ly a, 01 6300A, 
and 01 at the disk radii of 1AU, 10AU, and 

100AU with decreasing gas temperature. These 
properties are qualitatively the same even if we 
use the different dust models in §2.11 

Furthermore, we plot in Figure [7] the vertical 
profiles of the X-ray ionization rates, (x, defined 
in equation (fl"8j). at disk radii of 1AU, 10AU, and 
100AU, where the disk is irradiated by both of X- 
ray and UV radiation from the central star. The 
radial (Cx,i?,; dashed lines) and vertical (Cx,z! dot- 
ted lines) components, which are calculated by 
substituting Fx,r and Fx,z of equations ([9]) and 
([TO]) into equation (fT8|) and satisfy (x = (x,R + 
Cx,z, are also plotted for comparison. In addition, 
the ionization rates caused by interstellar cosmic- 
ray, £cr are plotted in dot-dashed lines, which are 
estimated as 



OR 



CcR,oexp[-S(z)/x 



CR 



(23) 



where we adopt Ccr.o = 1 x 10~ 17 s -1 and the 
attenuation coefficient of xcr = 96 g cm -2 (Ume- 
bayashi & Nakano 1981). The surface density 
is calculated as E(z) = p{z')dz' . The fig- 
ure shows that at the disk surface the ionization 
rates due to X-rays from the central star are much 
higher than those due to interstellar cosmic-rays, 
while near the disk midplane the former is much 
lower than the latter. This is because X-ray at- 
tenuation is larger than that of cosmic-rays and 
because the Compton scattering of X-ray radia- 
tion is inefficient (see Fig. [1] and Igea & Glassgold 
1999). 
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Fig. 6. — The vertical profiles of the cooling and 
heating rates at the disk radii of (a) x = 1 AU, (b) 
10 AU, and (c) 100 AU for the irradiation model of 
X-rays + UV and the dust model A with a max — 
lOfim. The X-ray or FUV heating dominates the 
heating process, while the radiative cooling (Ly 
a, OI 6300A, and OI 63^m for £ = 1, 10, and 
100 AU) and the dust-gas collision dominate the 
cooling process at the surface layer and near the 
midplane, respectively. 
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Fig. 7. — The vertical profiles of the X-ray ioniza- 
tion rate at the disk radii of x = 1, 10, and 100 
AU for the irradiation model of X-rays + UV and 
the dust model A with a max = lOjiim. The solid, 
dashed, and dotted lines show the total rate, radial 
and vertical components, respectively. The ioniza- 
tion rate due to interstellar cosmic-rays are plotted 
as a dot-dashed line for comparison. The ioniza- 
tion rates by the X-rays from the central star are 
much lower than those by the interstellar cosmic- 
ray near the disk midplane due to the inefficient 
Compton scattering of the X-ray radiation. 
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3.2.2. Effect of Dust Evolution 

In Figure [5] we plot the vertical gas tempera- 
ture profiles for various dust models at disk radii 
of x = 1, 10, and 100 AU for the case of a disk 
heated by both X-rays and UV radiation. The 
profiles for dust model A with different maximum 
dust radii of a max = 10/zm (solid lines), 1mm 
(dashed lines), and 10cm (dot-dashed lines) are 
plotted together in Figure [8^,. The profiles for 
dust model B at 10 6 yr after the calculation starts 
are plotted in Figure [5Jd for the quiescent (solid 
lines) and turbulent (dot-dashed lines) disks. The 
profiles calculated by using the dense cloud dust 
model (which is the initial condition of the calcu- 
lation for the dust evolution and the dust model 
A with a max = 10/mi) are also plotted together 
in thin solid lines for comparison. The thin dot- 
ted lines in the figures are the dust temperature 
profiles. Figure [5£i shows that as the dust parti- 
cles grow and the total surface area of dust grains 
(/dust) decreases (see §2.1.3|) . the gas temperature 
at the disk surface drops because the grain pho- 
toelectric heating rate decreases (e.g., Aikawa & 
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Fig. 8. — The vertical temperature profiles of dust 
(thin dotted lines) and gas at the disk radii of 
x = 1, 10 and 100AU for (a) dust model A with 
flmax = 10/im (solid lines), 1mm (dashed lines), 
and 10cm (dot-dashed lines), and (b) model B in 
quiescent (solid lines) and turbulent (dashed lines) 
disks. The profiles for the dense cloud dust model 
(dust model A with a max = 10/im) are plotted 
here as thin solid lines. The irradiation model of 
X-rays + UV is used here. As the dust particles 
grow or settle toward the disk midplane the gas 
temperature in the middle layer and in the outer 
disk decreases due to the decrease of the grain pho- 
toelectric heating rate. Meanwhile, the dust and 
gas temperatures near the midplane increases due 
to smaller grain opacity and greater penetration 
of the irradiation from the central star. 



15 



Nomura 2006). The gas temperatures in the in- 
ner disk (x ~ 1 AU) and in the surface layer at 
x ~ 10 AU do not change because the X-ray heat- 
ing dominates in these regions. The dust tempera- 
ture at the disk surface decreases slightly with dust 
growth. Figure [SJd shows that clear differences ap- 
pear in the gas and dust temperatures between 
the dust model of dense clouds (dust model A with 
Omax = 10/xm) and the models with the dust evolu- 
tion in both quiescent and turbulent disks. For the 
models with dust evolution the gas temperature in 
the middle layer and the outer region of the disks 
drops owing to the decrease of /dust (see §2.1.3l and 
§3.1|) . while the dust and gas temperatures near 
the midplane increase because the grain opacity 
decreases and the irradiation from the central star 
can penetrate deeper in to the disks. At x ~ 1AU 
heating via irradiation dominates even near the 
midplane for the models with the dust evolution, 
whereas the viscous heating is dominant for the 
dense cloud dust model. The differences between 
the quiescent and turbulent disks are small be- 
cause the profiles of /dust are similar, especially in 
the surface layer. The dust and gas temperatures 
very close to the midplane are slightly higher in 
the turbulent disk due to the higher collision rate 
between the dust particles, which results in lower 
/dust and grain opacity (see £ 13. lj) . 

Dust growth and settling are also expected to 
impact on the X-ray heating rates and the gas tem- 
perature profile through the change in the pho- 
toionization cross section, <7 p h, part of which is 
contributed by heavy elements in dust grains (e.g., 
Glassgold et al. 1997; Wilms et al. 2000). Here 
we check the effects by simply adopting an ex- 
treme case, specifically that the contribution of 
dust grains to the cross section is negligible for 
dust model A with a max = 10cm and for dust 
model B. We modify the cross section in Malony 
et al. (1996) by simply assuming that the contri- 
bution by heavy elements in gas phase is about 
60%, on average, of the total cross section (Wilms 
et al. 2000). This modification of the cross sec- 
tion makes the gas temperature higher or lower 
by a factor of 2 at most at the disk surface, where 
the X-ray heating is dominant. At large z the gas 
temperature becomes slightly lower due to the de- 
crease of the photoionization rate, while at smaller 
z, where the influence of attenuation is more im- 
portant, the temperature becomes a bit higher ow- 
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Fig. 9. — The vertical gas density profiles at the 
disk radii of x = 1, 10 and 100AU for dust model 
B in quiescent {solid lines) and turbulent {dashed 
lines) disks. The profiles for the dense cloud dust 
model are also plotted together in thin solid (x = 
1AU), dotted (x = 10AU), and dot-dashed (x 
1 Of ). A I ' lines. The disks are puffed out more for 
the models with the dust evolution due to higher 
gas temperatures at the disk midplane and higher 
disk scale height. 

ing to the decrease of the attenuation coefficient, 
which results in the relatively stronger X-ray ra- 
diation field (e.g., Glassgold et al. 1997). The gas 
density at the disk surface also becomes higher or 
lower by a factor of 3 at most, according to the 
change in the gas temperature. The variation of 
the photoionization cross section due to the dust 
evolution will also slightly affect the level popula- 
tions and the line emission of molecular hydrogen 
through the changes in the thermal collision and 
the X-ray pumping rates. When the FUV heating 
or the UV pumping process dominates, however, 
the changes will be small. In the following sections 
we neglect these effects for simplicity. 

As the dust particles evolve in the disk, the gas 
density profile also changes since it is related to the 
gas temperature profile. In Figure [9] we plot the 
gas density profiles in the vertical direction at the 
disk radii of x = 1, 10, and 100AU, which are cal- 
culated by using dust model B in quiescent {solid 
lines) and turbulent [dashed lines) disks. The 
profiles for the dense cloud dust model are also 
plotted together as thin solid (x = 1AU), dotted 
(x = 10AU), and dot-dashed (x = 100AU) lines 
for comparison. The gas densities in the models 
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Fig. 10. — The contour plots of the gas temperature (solid lines) and density (dotted lines) distributions in 
the z/x vs. z plane for dust model A with (a) a max = 10/im, (b) 1mm, and (c) 10cm, and model B in (d) 
quiescent and (e) turbulent disks. The irradiation model of X-rays + UV is used here. 



17 



with the dust evolution are lower at the disk mid- 
plane and higher at the disk surface than those 
for the dense cloud dust model because of higher 
gas temperatures at the midplane and higher disk 
scale height (H = c sQ /£Ik)- In Figure [TU] we 
present the contour plots of the resulting gas tem- 
perature (solid lines) and density (dashed lines) 
profiles in the z/x vs. x plane. The contour levels 
are taken as T = 30, 100, 300, 1000, and 3000K, 
and p = 10~ 16 , 10~ 14 , 10~ 12 , and 10" 10 g cm" 3 . 
Dust model A with (a) a max = 10/mi, (b) 1mm, 
and (c) 10cm, and model B in (d) quiescent and 
(e) turbulent disks are used in these figures. 

Making use of the obtained density and temper- 
ature distributions, we calculate the continuum ra- 
diation of thermal dust emission from the disks by 
solving the radiative transfer equation, simply as- 
suming that the disks are face on to an observer. 
The resulting infrared (IR) spectra basically re- 
produce the median spectral energy distribution 
(SED) observed toward classical T Tauri stars 
(CTTSs) in the Taurus- Auriga molecular cloud 
(D'Alessio et al. 1999, 2006) for the dust model 
A with a max = 10/mi. For the models with larger 
maximum dust radii, the resulting IR dust emis- 
sion is weaker than the median SED by a factor of 
about 4 at the most (e.g., D'Alessio et al. 2001). 
The thermal dust emission from the disks for dust 
model B in both quiescent and turbulent disks also 
reproduces the median observed SED if we adjust 
the inner disk radii. These models do not repro- 
duce the flux deficits relative to the median SED 
in the near-IR to the mid-IR wavelength bands 
which are observed toward several CTTSs, includ- 
ing TW Hya (e.g., Calvet et al. 2002; Bergin et al. 
2004). However, the disk structure beyond sev- 
eral tens of AU where molecular hydrogen lines 
are mainly emitted (see Paper I) will be almost 
unaffected even if we were to modify the structure 
of inner disk within several AU in order to force 
the thermal dust emission to account for the flux 
deficits, because only a limited region close to the 
midplane of the outer disk can be shadowed by the 
inner disk since the disk has a flared structure. 

Finally, in Figure fTTI we show another effect of 
the dust evolution due to the change in grain opac- 
ity. The figure shows the profiles of the integrated 
FUV radiation fields for the energy range of 6eV 
< hv < 13eV in the vertical direction at the disk 
radii of x = 1, 10, and 100 AU. The profiles for 
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Fig. 11. — The vertical profiles of the integrated 
FUV radiation fields (6eV < hv < 13eV) at the 
disk radii of x = 1, 10 and 100AU for dust models 
A (a) and B (b). The profiles are plotted in the 
same way as in Figure [8] The FUV radiation from 
the central star penetrates deeper in the disks as 
the dust particles evolve and the grain opacity de- 
creases. 

dust models A and B are plotted in Figure [TTk 
and b, respectively, in the same way as Figure [8] 
The figures show that as the dust particles evolve 
in the disk and /dust decreases, the FUV radia- 
tion from the central star penetrates deeper in the 
disk due to the decrease of grain opacity (see also 
Paper I for the FUV radiation fields in disks). 

3.3. Level Populations of Molecular Hy- 
drogen 

Making use of the physical properties of the 
disks obtained in the previous subsections, we cal- 
culate the level populations of molecular hydrogen 
in the disks by solving the equations for statistical 
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equilibrium f £|2.4p . The effects of X-ray irradia- 
tion and dust evolution on the level populations 
are discussed in the following. 

3.3.1. Effect of X-rays 

In Figure [12] we plot the resulting level popu- 
lations of molecular hydrogen for the models in 
which the disk is irradiated by (a) both X-ray and 
UV radiation from the central star, (b) UV ra- 
diation only, and (c) X-ray radiation only. Dust 
model A with maximum dust radius of a max = 
10/im (see §2.1|) is used in this sub-subsection. 
The filled diamonds show the column densities of 
molecular hydrogen in each ro- vibrational level as 
a function of the level energy. The column densi- 
ties are calculated by integrating the number den- 
sity of molecular hydrogen in each level along the 
vertical direction at a disk radius of 50AU. The 
level populations in local thermodynamic equilib- 
rium (LTE) are shown as dashed lines. The figure 
shows that if we take into account UV irradia- 
tion from the central star, the gas temperature 
becomes high enough for the collisional excitation 
process to be very efficient, and the level popula- 
tions in lower energy levels are in LTE distribution 
as a result. Meanwhile, if the disk is irradiated by 
X-rays only and the gas is cold enough, the pop- 
ulations are not in LTE due to the X-ray pump- 
ing process. It suggests that we may be able to 
observe molecular hydrogen transitions excited by 
X-ray pumping toward those protoplanetary disks 
whose central stars have strong X-ray and weak 
UV radiation. 

3.3.2. Effect of Dust Evolution 

We now discuss the effect of dust evolution on 
the level populations of molecular hydrogen in a 
disk irradiated by both X-rays and UV radiation. 
Figure [13] is the same as Figure [12] but for differ- 
ent dust models. Dust model A with maximum 
dust radii of a max = (a) 10/im, (b) 1mm, and 
(c) 10cm and (d) model B are used in these fig- 
ures. Figures IT3r-c show that as the dust particles 
grow, the level populations of molecular hydrogen 
change from LTE to non-LTE distributions. This 
is because with increasing dust size and decreas- 
ing /dust, the gas temperature drops due to the 
decrease of grain photoelectric heating rate, and 
the collisional excitation process becomes less effi- 
cient. In addition, the UV radiation from the 
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Fig. 12. — The level populations of molecular hy- 
drogen at a disk radius of 50 AU (filled diamonds) 
for the irradiation models of (a) X-rays + UV, (b) 
UV only, and (c) X-rays only. The LTE distribu- 
tions are plotted as dashed lines. Dust model A 
with a max = 10/im is used here. The populations 
are in LTE in lower energy levels when the disk is 
irradiated by strong UV radiation, while they are 
controlled by X-ray pumping if the UV irradiation 
is weak and X-ray irradiation is strong. 
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Fig. 13. — The level populations of molecular hydrogen at a disk radius of 30 AU for dust model A (filled 
diamonds) with o max = (a) 10/im, (b) 1mm, and (c) 10cm, and (d) model B in quiescent (asterisks) and 
turbulent (open squares) disks. The irradiation model of X-rays + UV is used here. The level populations 
change from LTE to non-LTE as dust particles grow or settle toward the disk midplane, and since the gas 
temperature drops while the UV photons in the disk increase. 
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central star can penetrate deeper in the disk due to 
the decrease of grain opacity, and the UV pumping 
process becomes more efficient. X-ray pumping is 
not the dominant process if we take into account 
the UV irradiation from the central star. In Fig- 
ures IT3t i the level populations with dust model 
B in quiescent (asterisks) and turbulent (open 
squares) disks are plotted together. The popula- 
tions in quiescent and turbulent disks arc almost 
identical because of similar physical properties of 
the disks (see £13 . 2 . 2|) . The figure shows that if we 
take into account dust evolution, the level popu- 
lations are in non-LTE distributions both in qui- 
escent and turbulent disks because of low /dust, 
which results in low gas temperature and high UV 
radiation fields, while they are in LTE distribu- 
tions with the dust model of dense clouds (dust 
model A with a max = 10/xm). 

3.4. Molecular Hydrogen Emission 

Making use of the physical properties of the 
disks and the level populations of molecular hydro- 
gen obtained in the previous subsections, we cal- 
culate the line emission from molecular hydrogen 
f ^2.4p . In the following we show the resulting line 
spectra in the near- and mid-infrared (NIR and 
MIR), and ultraviolet (UV) wavelength bands, 
and present line ratios, using various dust and ir- 
radiation models. 

3.4-1- Line Spectra 

Figures HU [15] and [16] show that the resulting 
line spectra in the NIR, MIR, and UV wavelength 
bands, respectively. Dust model A with a max = 
(a) 10/zm, (b) 1mm, and (c) 10cm, and (d) model 
B in a quiescent disk are used in these figures. The 
spectra from a turbulent disk are not plotted in 
these figures as they are almost identical to those 
in a quiescent disk. 

In Figure [14] and Table [T] (upper rows) we 
present the ro-vibrational line fluxes of molecu- 
lar hydrogen in the NIR wavelength band. These 
show that as the dust particles grow, the transi- 
tion lines from higher vibrational energy levels be- 
come relatively stronger in this wavelength band. 
This is because the level populations change from 
LTE to non-LTE distributions due to the decrease 
of the gas temperature and the increasing impor- 
tance of UV pumping, so that the populations in 



higher vibrational energy levels become relatively 
larger as we have seen in the previous section (see 
Fig. [T3"|) . The line fluxes decrease with the dust 
evolution because the area of high temperature re- 
gion in the disk shrinks. In Figure [T5l and Table Q] 
(lower rows) we present the pure rotational tran- 
sition lines in the MIR wavelength band. In this 
case the transition lines from lower energy levels 
become relatively stronger. This is because the 
level populations in the ground vibrational state 
are in LTE for all dust models, and the popu- 
lations in lower energy levels become relatively 
larger as the gas temperature decreases with in- 
creasing dust size. The line fluxes decrease with 
dust evolution for the same reason as the NIR 
lines. The MIR flux for dust model B, however, 
does not decrease so much because /dust is not 
very small in the outer disk (see Fig. [2^) . The line 
fluxes from lower energy levels are rather stronger 
than those in model A since the gas temperatures 
near the midplane, to which the MIR lines are sen- 
sitive, are higher for model B (see ij3.2.2[) . In Table 
[T]we also list the infrared line fluxes for different 
irradiation models, calculated by using dust model 
A with a max = 10/^m. For the irradiation model of 
X-rays only, the intensity of NIR lines is similar to 
that for the irradiation model of X-rays + UV and 
dust model A with a max = 1mm. The MIR lines 
are relatively weaker because the gas temperature 
at the outer disk is not so high if the irradiation 
source is X-rays only (see Fig. [5]). For the irra- 
diation model of UV only, the NIR and MIR line 
fluxes are a bit weaker than those for the irradia- 
tion model of X-rays + UV and the dust model A 
with a max = 10/im since the gas temperature for 
the latter model is higher due to the X-ray heating. 
The intensity of emission lines from lower energy 
levels in the MIR is similar between the former 
and the latter models because X-ray heating does 
not affect the gas temperature in the outer disk. 

Comparing our results for the 2.12 fim v = 
1 — ► <S*(1) transition to the observational data to- 
wards TW Hya of 1.0 x 10~ 15 ergs s^ 1 cm~ 2 (Bary 
et al. 2003), dust model A with a max = 1mm 
(/dust = 0.1) seems to be most suitable. We may 
need larger amount of small dust grains than that 
predicted in dust model B in order to reproduce 
the observed 2.12 fim line flux. The calculated 
fluxes of the MIR lines are consistent with the up- 
per limits of the ground based observations 
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Fig. 14. — The near-infrared < A < spectra of ro-vibrational transition lines of molecular 

hydrogen from the disks for dust model A with a max = (a) 10^m, (b) 1mm, and (c) 10cm, and (d) model B 
in a quiescent disk. The irradiation model of X-rays + UV is used here. The distance to the disk is set to 
be d = 56pc. The lines from higher energy levels become relatively stronger as the dust particles evolve and 
the level populations change from LTE to non-LTE. 
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Fig. 15. — Same as Fig. [21 but for the mid-infrared (5/im < A < 30/im) spectra of pure rotational emission 
lines. The lines from lower energy levels become relatively stronger as dust particles evolve and the gas 
temperature decreases. 
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Fig. 16. — Same as Fig. [HI but for the ultraviolet (1100A< A < 1800A) emission lines. The line fluxes of 
transitions originally pumped from higher energy levels seems to be relatively weaker as dust particles evolve 
and the gas temperature decreases. The lines are also affected by the strength of UV radiation field in the 
disk. 



24 



(Richter et al. 2002; Sako et al. 2005). 

In Figure [16] and Table [2] we present the line 
fluxes in the UV wavelength band. In the table 
we list the lines pumped by - 2 P(0), 0-2 P(l), 
1-2 P(5), 1-2 i?(6), 3-lP(14), and4-3P(5) 
transitions in the wavelength band of the strong 
Ly a emission from the central star (see e.g., Her- 
czeg et al. 2002; Paper I). The figure and table 
show that the fluxes of most lines decrease as dust 
particles evolve in the disk and the area of the 
high temperature region shrinks. The fluxes of 
some lines, however, do not decrease because the 
lines in the UV wavelength band, excited by UV 
photons, are affected very much by the strength of 
the UV radiation field in the disk, which becomes 
stronger with dust evolution (see Fig. [TTj) . The 
line fluxes for dust model B are relatively strong 
since the populations of molecular hydrogen in the 
energy levels in the ground electronic states, from 
which the lines are pumped, are larger. The pop- 
ulations in lower energy levels are large because 
the gas temperatures near the midplane are high, 
while those in higher energy levels, which are con- 
trolled by the UV pumping process, are large be- 
cause the UV irradiation from the central star pen- 
etrates deeper in the disk (see Figs. 151 and [TTj) . We 
note that in this work the fluxes of lines in the 
UV wavelength band caused by X-ray pumping 
are not calculated. Such pumping will affect the 
strength of weak emission lines of molecular hy- 
drogen (e.g., Bergin et al. 2004), but the effect on 
the flux of the strong lines, listed in Tabled which 
are pumped by the strong Ly a line emission, will 
be negligible. 

The calculated UV line fluxes are 10~ 15 ~ 
10~ 14 ergs s _1 cm~ 2 , which are consistent with 
the observations towards TW Hya (Herczeg et al. 
2002) to zeroth order but have discrepancies in 
details. This could be because the UV line fluxes 
depend not only on the density and temperature 
profiles of protoplanetary disks, but also on the 
strength and shape of the Ly a line irradiated from 
the central star, as shown by Herczeg et al. (2002, 
2004; see also Paper I). A simple single Gaussian 
profile is used for the Ly a line profile in this pa- 
per, while the actual line profile seems to be influ- 
enced by wind absorption (see Herczeg et al. 2002, 
2004). Thus, a more detailed analysis of the Ly a 
line profile, beyond the scope of this work, will be 
needed in order to fit the observed UV line fluxes 



in more detail. 

Our results show that molecular hydrogen emis- 
sion is strong and will be easier to observe toward 
those disks whose central stars have strong UV 
and X-ray radiation. In addition, if the disk con- 
tains relatively large amount of small dust grains, 
the volume of hot gas in the disk will be larger, 
and emission lines will be stronger. Therefore, the 
disks which have an observable signature of the 
presence of small dust grains, such as strong 10/im 
silicate emission, could be good targets in which 
to observe molecular hydrogen lines. 

3.4.2. v = 1 -> S(l)/v = 2 1 5(1) Line 
Ratio 

Finally, we discuss the effect of the dust evo- 
lution on a particular line ratio, the v = 1 — » 
S(l)/v = 2 — ► 1 5(1) ratio, which is often used 
as a probe of physical properties of astronomical 
objects. In Figure [17] the resulting line ratios for 
various irradiation models and dust models are 
plotted. The lines with diamonds in the left hand 
side of the figure show the ratios calculated us- 
ing dust model A with a max = 10/im, 1mm, and 
10cm (/dust = 1.0, 0.1, and 0.01) and irradiation of 
both X-rays and UV (solid lines), UV only (dashed 
lines), and X-rays only (dot- dashed lines). They 
show that if the disk is irradiated by the UV ra- 
diation from the central star, this ratio becomes 
larger as the dust particles grow and the total sur- 
face area of dust grains, /dust, decreases. This 
is because the level populations of molecular hy- 
drogen change from LTE to non-LTE distributions 
due to the increase in the grain photoelectric heat- 
ing rate and the decrease of grain opacity, as dis- 
cussed in ^3.3.21 This effect appears to be more 
efficient if the disk is irradiated only by UV and 
the heating source of the gas at the surface disk is 
grain photoelectric heating only. If the disk is ir- 
radiated by X-rays only, the ratio does not change 
with the maximum dust size. We note, however, 
that we have neglected the change in the X-ray 
photoionization cross section due to the dust evo- 
lution. If we take it into account, it affects the 
line ratio slightly with the error being about 25% 
for the model with a max = 10cm in the extreme 
case that the contribution to the cross sections of 
heavy elements in the dust is negligible (see ij3.2.2l 
Wilms et al. 2000). Figure Q7| shows that the line 
ratio in the model with X-ray irradiation only is 
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slightly larger than that in the model with UV ir- 
radiation, although the gas temperature is lower. 
This occurs because the level populations are not 
in LTE, but affected by the X-ray pumping process 
(see M3.2. 1 [) . We also plot in the right hand side of 
Figurc[l7]thc resulting line ratios for dust model B 
in quiescent and turbulent disks in filled and open 
triangles, respectively. The ratio for the model 
with dust grains typical of dense cloud (model A 
with a max = 10/Ltm) is plotted as a filled circle. 
The results suggest that if the dust grains coag- 
ulate and settle towards the disk midplane, the 
ratio becomes substantially larger than the case 
in which dust grains do not evolve. So, our re- 
sults suggest that dust evolution in protoplanetary 
disks could be observable through this particular 
line ratio. Itoh et al. (2003) derived an upper 
limit to the line ratio of 0.26 from an observation 
toward a classical T Tauri star, LkHa 264, and 
Bary et al. (2007, in preparation) find an upper 
limit of ~0.2 toward a Herbig Be star, HD 97048; 
all the models we used in this paper almost satisfy 
these observational upper limits. 

4. Summary 

We have made a detailed model of physical 
structure of protoplanetary disks and calculated 
the level populations and line emission of molec- 
ular hydrogen, taking into account X-ray irradia- 
tion from the central star as well as dust growth 
and settling towards the disk midplane. We have 
followed the time evolution of the spatial and size 
distributions of dust particles in the disks by nu- 
merically solving the coagulation equation for set- 
tling dust particles. The resulting mass and total 
surface area of dust grains per unit gas volume is 
much smaller, except at the disk midplane, than 
those for the model in which the dense cloud dust 
grains are well mixed with the gas. At the disk 
surface the dust density normalized by the ini- 
tial value, Pdust/pdust,0j and the parameter for the 
total surface area of dust grains, /dust, are small 
due to the dust settling toward the disk midplane. 
Near the disk midplane the parameter /dust be- 
comes much smaller since small dust particles are 
removed by dust coagulation. 

We have studied the effects of X-ray irradiation 
on the physical structure of the disks and found 
that the X-ray irradiation is the dominant heating 
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Fig. 17.— The v = 1 -> S(l)/v = 2 1 S(l) 
line ratios for the irradiation models of X-rays + 
UV (solid lines), UV only (dashed lines), and X- 
rays only (dot- dashed lines) and dust model A with 
flmax = 10/im, 1mm, and 10cm (/dust = 1-0, 0.1, 
and 0.01) (in the left-hand panel). The ratios for 
the irradiation model by X-rays + UV with dust 
model B in quiescent (open triangle) and turbu- 
lent (filled triangle) disks and the dense clouds 
(dust model A with a max = 10/im; filled circle) 
are plotted in the right-hand panel. If the disk is 
irradiated by strong UV radiation, the line ratio 
increases with dust evolution as a consequence of 
the decrease in gas temperature and the increase 
in the UV radiation field strength . 

source in the inner region and in the surface layers 
of the disk. FUV heating dominates in the middle 
layers and in the outer region of the disk. This is 
because the FUV radiation is scattered efficiently 
by dust grains, while the Compton scattering of X- 
ray radiation is inefficient in the energy range of 
E < lkcV in which T Tauri stars mainly emit X- 
ray radiation. We found that the ionization rate 
caused by X-rays is much smaller than that due 
to interstellar cosmic-rays near the disk midplane 
because of the relatively large attenuation and in- 
efficient scattering of X-rays. 

The dust evolution in the disks affects the phys- 
ical disk structure, especially the gas temperature 
at the disk surface and the FUV radiation field 
within the disk. As the dust particles grow or set- 
tle towards the disk midplane, the gas temperature 
in the middle layers and the outer disk decreases 
because the grain photoelectric heating which is 
induced by FUV radiation becomes less efficient. 
Meanwhile, the FUV radiation from the central 
star penetrates deeper into the disk due to the de- 
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crease of grain opacity. 

Furthermore, making use of the obtained phys- 
ical structure of the disks, we calculated the level 
populations of molecular hydrogen in the ground 
electronic state. Our results show that if the cen- 
tral star has strong X-ray and weak UV radia- 
tion, the level populations arc controlled by X- 
ray pumping. Otherwise, the level populations 
are mainly controlled by thermal collisions or UV 
pumping, depending on the dust properties in the 
disk. As the dust particles evolve in the disk, the 
level populations change from LTE to non-LTE 
distributions since collisional excitation becomes 
less efficient due to the decrease of the gas tem- 
perature at the disk surface while UV pumping 
becomes more efficient owing to the stronger UV 
radiation field in the disk. 

Finally, using these level populations, we cal- 
culated the line emission of molecular hydrogen 
from the disk. The ro-vibrational line spectra 
in the near-infrared wavelength band show that 
the emission lines from high energy levels become 
relatively stronger as the dust particles evolve in 
the disk. Again, this is due to level populations 
changing from LTE to non-LTE and the popula- 
tions in higher vibrational energy levels becom- 
ing relatively larger. For the pure rotational line 
spectra in the mid-infrared wavelength band, it is 
the emission lines from lower energy levels which 
become relatively stronger. This is because the 
level populations in the ground vibrational state 
are in LTE and the populations in lower energy 
levels become relatively larger with dust evolution 
and decreasing dust surface area which results in 
lower gas temperatures. For transitions in the UV 
wavelength band, the dependence on the dust evo- 
lution is not so straightforward because the line 
fluxes decrease as the area of high temperature 
region shrinks, while they increase as UV irradia- 
tion from the central star penetrates deeper in the 
disks. Basically, the line fluxes which originate 
from pumping from higher energy levels seem to 
be relatively weaker as the dust particles evolve 
and the gas temperature decreases. Our results 
suggest that infrared line ratios of molecular hy- 
drogen could be a useful probe of dust evolution in 
protoplanetary disks. If the dust particles evolve, 
the v = 1 -> S(l)/v = 2 -> 1 5(1) line ra- 
tio, for example, becomes clearly larger than that 
for the dense cloud dust model (without the dust 



evolution). Further observations of the line ratios 
of molecular hydrogen could provide some con- 
straints on the dust evolution model in protoplan- 
etary disks. 

We are grateful to an anonymous referee for 
his comments which improve the clarity of our 
discussion. We would like to thank J.S. Bary, 
D.A. Weintraub, M.J. Richter, and M.A. Bitner 
for giving us information about unpublished obser- 
vational data, Y. Itoh and T. Takeuchi for fruitful 
comments, and T. Matsuda for useful help for nu- 
merical calculations. This work is supported by 
"The 21st Century COE Program of Origin and 
Evolution of Planetary Systems" and the Grants- 
in-Aid for Scientific Research 17039008, 17540217, 
and 18026006 in MEXT. Astrophysics at Queen's 
University Belfast is supported by PPARC. H. N. 
and M. T. acknowledge financial supported from 
the Japan Society for the Promotion of Science. 

REFERENCES 

Adachi, I., Hayashi, C, & Nakazawa, K. 1976, 
Prog. Thcor. Phys., 56, 1756 

Aikawa, Y. & Herbst, E. 1999, A&A, 351, 233 

Aikawa, Y. & Herbst, E. 2001, A&A, 371, 1107 

Aikawa, Y. & Nomura, H. 2006, ApJ, 642, 1152 

Ardila, D.R., Basri, G., Walter, F.M., Valenti, 
J.A., & Johns-Krull, CM. 2002, ApJ, 566, 1100 

Bary, J.S., Weintraub, D.A., & Kastner, J.H. 

2002, ApJ, 576, L73 

Bary, J.S., Weintraub, D.A., & Kastner, J.H. 

2003, ApJ, 586, 1136 

Beckwith, S., Gatley, I., & Pcrsson, S.E. 1978, 
ApJ, 219, L33 

Beckwith, S.V.W., & Sargent, A.I. 1993, in Proto- 
stars and Planets III, ed. E.H. Levy & J. Lunine 
(Tuscon, Univ. Arizona Press), 521 

Bergin, E.A., Aikawa, Y., Blake, G.A., & van 
Dishoeck, E.F. 2007, in Protostars and Plan- 
ets V, cd. B. Reipurth, D. Jewitt, and K. Keil 
(Tuscon: Univ. Arizona Press), in press 

Bergin, E.A., Calvet, N., D'Alessio, P. & Herczeg, 
G. J. 2003, ApJ, 591, L159 



27 



Bergin, E.A. et al. 2004, ApJ, 614, L133 

Black, J. H., & van Dishoeck, E. F. 1987, ApJ, 
322, 412 

Bohrcn, C. F., & Huffman, D. R. 1983, Absorption 
and scattering of light by small particles (New 
York: Wiley) 

Brown, A., Millar, T.J., Williams, P.M., & Zealey, 
W.J. 1983, MNRAS, 203, 785 

Burton, M. G., Bulmcr, M., Moorhouse, A., Ga- 
balle, T. R., & Brand, P. W. J. L. 1992, MN- 
RAS, 257, 1 

Calvet, N., D'Alessio, P., Hartmann, L., Wilncr, 
D., Walsh, A., & Sitko, M. 2002, ApJ, 568, 1008 

Carr, J.S. 1989, ApJ, 345, 522 

Cazaux, S. & Ticlcns, A.G.G.M. 2002, ApJ, 575, 
L29 

Cazaux, S. & Tielens, A.G.G.M. 2004, ApJ, 604, 
222 

Costa, V. M., Lago, M. T. V. T., Norci, L., & 
Meurs, E. J. A. 2000, A&A, 354, 621 

Cuzzi, J. N., Dobrovolskis, A. R., & Champney, J. 
M. 1993, Icarus, 106, 102 

D'Alessio, P., Calvet, N., & Hartmann, L. 2001, 
ApJ, 553, 321 

D'Alessio, P., Calvet, N., Hartmann, L., Franco- 
Hernandez, R., Scrvin, H. 2006, ApJ, 638, 314 

D'Alessio, P., Calvet, N., Hartmann, L., Lizano, 
S., & Canto, J. 1999, ApJ, 527, 893 

Duchcnc, G., McCabe, C, Ghcz, A. M., & Mac- 
intosh, B. A. 2004, ApJ, 606, 969 

Dullemond, C. P. & Dominik, C. 2004, A&A, 421, 
1075 

Dullemond, C. P. & Dominik, C. 2005, A&A, 434, 
971 

Dullemond, CP., Hollcnbach, D., Kamp, I., & 
D'Alessio, P. 2007, in Protostars and Planets V, 
ed. B. Rcipurth, D. Jewitt, and K. Keil (Tus- 
con: Univ. Arizona Press), in press 



Dullemond, C. P. & Turolla, R. 2000, A&A, 360, 
1187 

Dutrcy, A., Guilloteau, S., & Guelin, M. 1997, 
A&A, 317, L55 

Dutrcy, A., Guilloteau, S., & Ho, P. 2007, in 
Protostars and Planets V, ed. B. Rcipurth, D. 
Jewitt, and K. Keil (Tuscon: Univ. Arizona 
Press), in press 

Feigelson, E. D., & Montmerlc, T. 1999, ARA&A, 
37, 363 

Getman, K. V., et al. 2005, ApJS, 160, 319 

Glassgold, A. E., Najita, J., & Igea, J. 1997, ApJ, 
480, 344 

Glassgold, A.E., Najita, J., & Igea, J. 2004, ApJ, 
615, 972 

Glassgold, A.E. & Najita, J.R. 2001, in ASP Conf. 
Ser. 244, Young Stars near Earth: Progress and 
Prospects, ed. R. Jayawardhana & T. P. Greene 
(San Francisco: ASP), 251 

Gorti, U, & Hollcnbach, D. 2004, ApJ, 613, 424 

Gredcl, R., & Dalgarno, A. 1995, ApJ, 446, 852 

Gredel, R., Lepp, S., Dalgarno, A. & Herbst, E. 
1989, ApJ, 347, 289 

Hasegawa, T., Gatley, I., Garden, R. P., Brand, P. 
W. J. L., Ohishi, M., Hayashi, M., & Kaifu, N. 
1987, ApJ, 318, L77 

Hcrbig, G.H. & Goodrich, R.W. 1986, ApJ, 309, 
294 

Herbst, W., Herbst, D.K., Grossman, E.J., & We- 
instein, D. 1994, AJ, 108, 1906 

Herczeg, G.J., Linsky, J.L., Valcnti, J. A., Johns- 
Krull, CM., & Wood, B.E. 2002, ApJ, 572, 310 

Herczeg, G.J., Linsky, J.L., Walter, F.M., Gahm, 
G.F., & Johns-Krull, CM. 2006, ApJS, 165, 
256 

Herczeg, G.J., Wood, B.E., Linsky, J.L., Valenti, 
J.A., & Johns-Krull, CM. 2004, ApJ, 607, 369 

Hollenbach, D. & McKee, C.F. 1989, ApJ, 342, 
306 



28 



Igea, J. & Glassgold, A. E. 1999, ApJ, 518, 848 

Imanishi, K., Nakajima, H., Tsujimoto, M., 
Koyama, K., & Tsuboi, Y. 2003, PASJ, 55, 653 

Itoh, Y. et al. 2003, PASJ, 54, 963 

Itoh, Y., Sugitani, K., Ogura, K., & Tamura, M. 
2003, PASJ, 55, L77 

Jonkhcid, B., Faas, F. G. A., van Zadclhoff, G.-J., 
& van Dishoeck, E. F. 2004, A&A, 428, 511 

Jonkheid, B., Kamp, I., Augereau, J.-C, & van 
Dishoeck, E. F. 2006, A&A, 453, 163 

Jonkhcid, B., Dullcmond, CP., Hogerheijde, 
M.R., fe van Dishoeck, E. F. 2007, A&A, in 
press ( |astro-p h/0611223) 

Kaastra, J. S. 1992, An X-Ray Spectral Code for 
Optically Thin Plasmas (Internal SRON-Leidcn 
Report, updated version 2.0) 

Kamp, I. & Dullemond, C. P. 2004, ApJ, 615, 991 

Kamp, I., Dullemond, CP., Hogerheijde, M., & 
Enriqucz, J.E. 2005, in IAU Symp. 231, Astro- 
chemistry: Recent Successes and Current Chal- 
lenges, ed. D. C. Lis, G. A. Blake, & E. Herbst 
(Cambridge: Cambridge Univ. Press), 377 

Kastner, J. H., Huencmocrder, D. P., Schulz, N. 
S., & Canizares, C. R. 2002, ApJ, 567, 434 

Kenyon, S. J., & Hartmann, L. 1987, ApJ, 323, 
714 

Kenyon, S. J., & Hartmann, L. 1995, ApJS, 101, 
117 

Kitamura, Y., Momose, M., Yokogawa, S., 
Kawabe, R., Tamura, M., & Ida, S. 2002, ApJ, 
581, 357 

Koyama, K., Maeda, Y., Ozaki, M., Ueno, S., Ka- 
mata, Y., Tawara, Y., Skinner, S., & Yamauchi, 
S. 1994, PASJ, 46, L125 

Kozasa, T., Blum, J., & Mukai, T. 1992, A&A, 
263, 423 

Landau, L.D., Akhiezer, A.I., & Lifshitz, E.M. 
1967, General Physics (Oxford; New York: 
Pergamon Press) 



Liedahl, D. A., Osterheld, A. L., & Goldstein, W. 
H. 1995, ApJ, 438, L115 

Maloney, P.R., Hollenbach, D.J., & Tielens, 
A.G.G.M. 1996, ApJ, 466, 561 

Markwick, A. J., Ilgner, M., Millar, T. J., & Hen- 
ning, T. 2002, A&A, 385, 632 

McMaster, W.H., Kerr Del Grande, N., Mallett, 
J.H., & Hubbell, J.H. 1969, in Compilation of 
X-Ray Cross Sections, Lawrence Livermore Na- 
tional Laboratory Report UCRL-50174, Sec II, 
Rev.l 

Mewe, R., Gronenschild, E. H. B. M., & van den 
Oord, G. H. J. 1985, A&AS, 62, 197 

Miyake, K., & Nakagawa, Y. 1993, Icarus, 106, 20 

Miyake, K., & Nakagawa, Y. 1995, ApJ, 441, 361 

Mizuno, H. 1989, Icarus, 80, 189 

Mizuno, H., Markiewicz, W. J., & V61k, H. J. 
1988, A&A, 195, 183 

Najita, J.R., Carr, J.S., Glassgold, A.E., & 
Valenti, J. A. 2007, in Protostars and Planets V, 
ed. B. Reipurth, D. Jewitt, and K. Keil (Tus- 
con: Univ. Arizona Press), in press 

Nakagawa, Y., Nakazawa, K., & Hayashi, C. 1981, 
Icarus, 45, 517 

Nakagawa, Y., Sekiya, M., & Hayashi, C. 1986, 
Icarus, 67, 375 

Nomura, H., 2002, ApJ, 567, 587 

Nomura, H. & Millar, T.J. 2005, A&A, 438, 923 
(Paper I) 

Nomura, H. & Nakagawa, Y. 2006, ApJ, 640, 1099 

Ossenkopf, V. 1993, A&A, 280, 617 

Ormel, C W., Spaans, M., & Tielens, A.G.G.M. 
2007, A&A, 461, 215 

Pirronello, V., Liu, C, Roser, J.E. & Vidali, G. 
1999, A&A, 344, 681 

Preibisch, T., et al. 2005, ApJS, 160, 401 

Qi et al. 2006, ApJ, 636, L157 



29 



Rettig, T., Brittain, S., Simon, T., Gibb, E., Bal- 
sara, D.S., Tilley, D.A., Kulesa, C. 2006, ApJ, 
646, 342 

Richter, M.J., Jaffe, D.T., Blake, G.A., & Lacy, 
J.H. 2002, ApJ, 572, L161 

Richter, M.J., Lacy, J.H., Jaffe, D.T., Greathouse, 
T.K., & Blake, G.A. 2007, ApJ, submitted 

Roddier, C, Roddier, F., Northcott, M. J., 
Graves, J. E., & Jim, K. 1996, ApJ, 463, 326 

Sako, S., Yamashita, T., Kataza, H., Miyata, T., 
Okamoto, Y.K., Honda, M., Fujiyoshi, T., & 
Onaka, T. 2005, ApJ, 620, 347 

Schmitt, W., Henning, T., & Mucha, R. 1997, 
A&A, 325, 569 

Shemansky, D. E., Hall, D. T., & Ajello, J. M. 
1985, ApJ, 296, 765 

Sheret, I., Ramsay Howat, S.K., & Dent, W.R.F. 
2003, MNRAS, 343, L65 

Stelzer, B. & Schmitt, J.H.M.M. 2004, A&A, 418, 
687 

Sternberg, A., & Dalgarno, A. 1989, ApJ, 338, 197 

Storey, P.J. & Hummer, D.G. 1995, MNRAS, 272, 
41 

Suttncr, G., & Yorkc, H. W. 2001, ApJ, 551, 461 

Takahashi, J. & Uehara, H. 2001, ApJ, 561, 843 

Takeuchi, T. & Lin, D.N.C. 2005, ApJ, 623, 482 

Tanaka, H., Himcno, Y., & Ida, S. 2005, ApJ, 625, 
414 

Tanaka, M., Hasegawa, T., Hayashi, S.S., Brand, 
P.W.J.L., Gatley, I. 1989, ApJ, 336, 207 

Thi, W.F. et al. 2001a, Nature, 409, 60 

Thi, W.F. et al. 2001b, ApJ, 561, 1074 

Thi, W.-F., van Zadelhoff, G.-J., & van Dishoeck, 
E. F. 2004, A&A, 425, 955 

Thi, W., van Dishoeck, E. F., Blake, G. A., van 
Zadelhoff, G., & Hogerheijde, M. R. 1999, ApJ, 
521, L63 



Tine, S., Lepp, S., Gredel, R., & Dalgarno, A. 
1997, ApJ, 481, 282 

Tsujimoto, M., Koyama, K., Tsuboi, Y., Goto, 
M., & Kobayashi, N. 2002, ApJ, 566, 974 

Umebayashi, T. & Nakano, T. 1981, PAS J, 33, 617 

Valenti, J. A., Johns-Krull, CM., & Linsky, J.L. 
2000, ApJS, 129, 399 

Verner, D. A. & Yakovlev, D. G. 1995, A&AS, 109, 
125 

Wagenblast, R. & Hartquist, T. W. 1988, MN- 
RAS, 230, 363 

Weidenschilling, S. J. 1977, MNRAS, 180, 57 

Weidenschilling, S. J. 1980, Icarus, 44, 172 

Weidenschilling, S. J. 1997, Icarus, 127, 290 

Weingartner, J. C, & Drainc, B. T. 2001, ApJ, 
548, 296 

Weintraub, D. A., Bary, J. S., Kastner, J. H., 
Shukla, S. J., & Chynoweth, K. 2005, in Pro- 
tostars and Planets V, ed. V. Mannings, et al. 
(LPI Contribution 1286; Houson: LPI), 8233 

Weintraub, D.A., Kastner, J.H., & Bary, J.S. 
2000, ApJ, 541, 767 

Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 
542, 914 

Zccho, T., Gurtler, A., Sha, X., Jackson, B. & 
Kiippers, J. 2002, J. Chem. Phys., 117, 8486 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.2. 



30 



Table 1 

Calculated infrared line fluxes of molecular hydrogen. 
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1.16 


2-0 S 


1) 


0.068 


0.182 


0.439 


0.051 





016 





020 


0.019 


1.19 


2-0 S 


0) 


0.013 


0.039 


0.091 


0.013 





005 





006 


0.006 


1.24 


2-0 Q 


(!) 


0.043 


0.143 


0.310 


0.056 





057 





067 


0.061 


1.24 


2-0 Q 


(2) 


0.014 


0.041 


0.095 


0.014 





008 





009 


0.009 


1.25 


2-0 Q 


(3) 


0.043 


0.116 


0.280 


0.033 





010 





014 


0.013 


1.83 


1-0 S 


5) 


0.587 


2.286 


4.849 


0.446 





029 





050 


0.049 


1.89 


1-0 S 


4) 


0.314 


1.555 


3.007 


0.287 





017 





027 


0.026 


1.96 


1-0 S 


3) 


1.326 


8.096 


14.475 


1.472 





079 





119 


0.111 


2.03 


1-0 S 


2) 


0.533 


3.857 


6.477 


0.719 





041 





058 


0.051 


2.12 


1-0 S 


l) a 


1.565 


12.968 


20.725 


2.526 





142 





169 


0.135 


2.22 


1-0 S 


0) 


0.357 


3.256 


5.023 


0.667 





040 





051 


0.037 


2.25 


2-1 S 


1) 


0.090 


0.242 


0.583 


0.068 





022 





028 


0.026 


2.35 


2-1 S 


0) 


0.019 


0.056 


0.130 


0.019 





011 





013 


0.012 


2.41 


1-0 Q 


(!) 


1.270 


12.299 


18.528 


2.707 





355 





484 


0.373 


2.41 


1-0 Q 


(2) 


0.395 


3.604 


5.560 


0.748 





066 





100 


0.077 


2.42 


1-0 Q 


(3) 


1.101 


9.128 


14.588 


1.780 





117 





185 


0.145 


2.44 


1-0 Q 


(4) 


0.298 


2.156 


3.621 


0.402 





023 





034 


0.029 


2.45 


1-0 Q 


5) 


0.648 


3.954 


7.069 


0.719 





039 





059 


0.055 


2.47 


1-0 Q 


(6) 


0.142 


0.702 


1.357 


0.129 





008 





012 


0.012 


2.50 


1-0 Q 


(7) 


0.260 


1.011 


2.145 


0.197 





013 





022 


0.021 


2.55 


2-1 Q 


(!) 


0.066 


0.223 


0.483 


0.088 





090 





110 


0.097 


2.56 


2-1 Q 


(2) 


0.022 


0.063 


0.147 


0.021 





012 





014 


0.013 


2.57 


2-1 Q 


(3) 


0.066 


0.176 


0.424 


0.050 





016 





021 


0.019 


6.11 


0-0 S 


6) 


0.368 


3.949 


5.791 


0.755 


0.050 


0.084 


0.058 


6.91 


0-0 S 


5) 


1.203 


17.346 


23.219 


3.928 


0.441 


0.880 


0.564 


8.02 


0-0 S 


4) 


0.375 


7.023 


8.753 


2.038 


0.427 


1.072 


0.731 


9.66 


0-0 S 


3) 


0.849 


19.835 


23.286 


8.035 


2.943 


7.160 


6.140 


12.27 


0-0 S 


2) b 


0.163 


4.271 


4.719 


2.751 


1 


688 


4.803 


4.711 


17.02 


0-0 S 


!)b,c 


0.209 


4.466 


4.606 


5.806 


5.148 


9.690 


9.900 


28.20 


0-0 S 


0) c 


0.020 


0.145 


0.142 


0.365 


0.532 


1.521 


1.717 



a Ground-bascd observations by Bary ct al. (2003) and Itoh ct al. (2003) give fluxes of (1.0 — 
15) X 10 -15 crg/s/cm 2 . 

Ground-based observations by Richter et al. (2002) and Sako et al. (2005) give fluxes of 
< 30 X 10~ 15 crg/s/cm 2 for the 12/im S(2) line and < 39 X 10~ 15 crg/s/cm 2 for the 17/xm S(l) 
line. 

c ISO observations by Thi ct al. (2001) give fluxes of (28-81) X 10~ 15 crg/s/cm 2 for the 17/xm 
S(l) line and (25 - 57) X lO" 15 erg/s/cm 2 for the 28Mm S(0) line. 
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Table 2 

Calculated and observed ultraviolet line fluxes of molecular hydrogen. 













(xl0~ 15 


erg/s/cm 2 


) 










Fl ino 


(xlO 15 erg/s/cm 


2 ) 




X 


Transition 


Model A (o m 




Model B 


Obs. a 


A 


Transition 


Model A (a m8 


x) 


Model B 


Obs. a 


(A) 






10/im 


1mm 


10cm 


Quics. 


Turb. 




(A) 




10/^m 


1mm 


10cm 


Quics. 


Turb. 




Pumped by 0-2 R(0) 


Pumped by 0-2 R(l) 


1161.7 


0-1 


R(0) 


0.436 


0.427 


0.291 


0.311 


0.300 




1162.2 


0-1 R(l) 


0.501 


0.512 


0.412 


0.428 


0.424 




1166.3 


0-1 


P(2) 


0.688 


0.700 


0.608 


0.638 


0.621 




1169.8 


0-1 P(3) 


0.681 


0.721 


0.637 


0.642 


0.633 




1217.3 


0-2 


R(0) 


1.728 


1.830 


1.245 


1.505 


1.585 




1217.7 


0-2 R(l) 


1.759 


1.909 


1.658 


1.786 


1.806 


9.1 


1222.0 


0-2 


P(2) 


2.719 


3.039 


2.623 


3.121 


3.398 




1225.6 


0-2 P(3) 


2.395 


2.683 


2.538 


2.715 


2.789 




1274.6 


0-3 


R(0) 


4.550 


4.478 


2.608 


3.630 


3.954 


27.4 


1275.0 


0-3 R(l) 


4.304 


4.918 


3.941 


4.506 


4.624 


24.6 


1279.6 


0-3 


P(2) 


7.690 


8.210 


5.606 


7.717 


8.722 


39.2 


1283.2 


0-3 P(3) 


5.950 


6.996 


6.092 


7.000 


7.443 


28.0 


1333.6 


0-4 


R(0) 


6.694 


5.909 


3.344 


5.087 


5.683 


42.8 


1333.9 


0-4 R(l) 


7.005 


7.487 


5.268 


6.252 


6.606 


7.9 


1338.7 


0-4 


P(2) 


12.495 


11.472 


7.185 


11.063 


12.584 


73.1 


1342.4 


0-4 P(3) 


10.007 


10.958 


8.251 


10.080 


11.168 


64.9 


1393.9 


0-5 


R(0) 


5.516 


4.866 


2.757 


4.199 


4.693 


35.3 


1394.1 


0-5 R(l) 


5.790 


6.177 


4.344 


5.159 


5.456 


52.4 


1399.1 


0-5 


P(2) 


10.308 


9.453 


5.923 


9.136 


10.389 


73.8 


1402.8 


0-5 P(3) 


8.277 


9.047 


6.807 


8.328 


9.235 


73.1 


1455.0 


0-6 


R(0) 


3.232 


2.819 


1.633 


2.563 


2.873 


20.8 


1455.2 


0-6 R(l) 


3.574 


3.686 


2.570 


3.111 


3.349 


30.9 


1460.4 


0-6 


P(2) 


6.141 


5.507 


3.483 


5.576 


6.290 


41.6 


1464.0 


0-6 P(3) 


5.148 


5.420 


4.030 


5.091 


5.744 


42.1 


1516.4 


0-7 


R(0) 


1.273 


1.097 


0.666 


1.107 


1.230 


21.2 


1516.5 


0-7 R(l) 


1.514 


1.491 


1.044 


1.325 


1.471 


21.2 


1521.8 


0-7 


P(2) 


2.467 


2.155 


1.383 


2.355 


2.597 


16.2 


1525.4 


0-7 P(3) 


2.202 


2.197 


1.617 


2.194 


2.515 


17.9 


1577.3 


0-8 


R(0) 


0.302 


0.258 


0.164 


0.291 


0.315 




1577.1 


0-8 R(l) 


0.378 


0.361 


0.259 


0.362 


0.412 




1582.6 


0-8 


P(2) 


0.592 


0.509 


0.330 


0.593 


0.637 




1586.0 


0-8 P(3) 


0.553 


0.531 


0.389 


0.590 


0.662 




1636.4 


0-9 


R(0) 


0.034 


0.029 


0.019 


0.035 


0.037 




1636.1 


0-9 R(l) 


0.044 


0.041 


0.030 


0.048 


0.053 




1641.7 


0-9 


P(2) 


0.067 


0.058 


0.038 


0.069 


0.073 




1644.8 


0-9 P(3) 


0.064 


0.061 


0.044 


0.073 


0.079 




1692.5 


0-10 


R(0) 


0.001 


0.001 


0.000 


0.001 


0.001 




1692.0 


0-10 R(l) 


0.001 


0.001 


0.001 


0.001 


0.001 




1697.6 


0-1C 


P(2) 


0.002 


0.001 


0.001 


0.002 


0.002 




1700.3 


0-10 P(3) 


0.002 


0.001 


0.001 


0.002 


0.002 










Pum 


pod by 1-2 P(5) 








Pumped by 1-2 R(6) 


1148.7 


1-1 


R(3) 


3.247 


3.386 


2.774 


2.701 


2.688 


4.6 


1162.0 


1-1 R(6) 


3.417 


3.415 


2.505 


2.559 


2.525 




1161.9 


1-1 


P(5) 


4.165 


4.369 


3.660 


3.631 


3.611 


10.9 


1183.4 


1-1 P(8) 


4.524 


4.591 


3.294 


3.467 


3.404 




1202.5 


1-2 


R(3) 


8.145 


8.716 


7.575 


7.909 


7.881 


11.3 


1215.8 


1-2 R(6) 


8.795 


9.284 


6.379 


6.986 


6.787 




1216.2 


1-2 


P(5) 


10.449 


11.295 


9.889 


10.442 


10.425 




1238.0 


1-2 P(8) 


11.411 


11.984 


7.394 


8.460 


8.162 


11.5 


1257.9 


1-3 


R(3) 


10.991 


11.831 


8.676 


9.699 


9.640 


18.1 


1271.2 


1-3 R(6) 


11.569 


11.074 


5.822 


7.064 


6.754 


14.1 


1272.0 


1-3 


P(5) 


13.849 


14.827 


10.630 


12.158 


12.127 


20.5 


1294.0 


1-3 P(8) 


13.476 


12.341 


6.279 


7.831 


7.423 


13.0 


1314.8 


1-4 


R(3) 


2.933 


2.700 


1.686 


2.061 


2.048 


12.2 


1327.8 


1-4 R(6) 


2.581 


2.181 


1.062 


1.362 


1.270 


6.1 


1329.3 


1-4 


P(5) 


3.536 


3.240 


2.022 


2.530 


2.489 


7.5 


1351.2 


1-4 P(8) 


2.689 


2.303 


1.129 


1.437 


1.347 


2.8 


1372.7 


1-5 


R(3) 


2.063 


1.892 


1.180 


1.449 


1.438 


3.2 


1385.2 


1-5 R(6) 


1.813 


1.530 


0.745 


0.956 


0.891 




1387.5 


1-5 


P(5) 


2.486 


2.272 


1.416 


1.776 


1.745 


7.1 


1409.2 


1-5 P(8) 


1.899 


1.621 


0.793 


1.011 


0.947 


2.2 


1431.2 


1-6 


R(3) 


12.444 


12.615 


8.287 


9.576 


9.548 


29.0 


1443.1 


1-6 R(6) 


12.251 


10.827 


5.388 


6.837 


6.427 


11.3 


1446.3 


1-6 


P(5) 


15.518 


15.441 


10.068 


12.060 


12.049 


44.2 


1467.4 


1-6 P(8) 


11.488 


10.895 


5.671 


7.012 


6.672 


17.6 


1489.8 


1-7 


R(3) 


15.769 


16.791 


11.439 


13.034 


12.946 


48.2 


1500.7 


1-7 R(6) 


16.461 


15.004 


7.589 


9.577 


9.034 


19.7 


1505.0 


1-7 


P(5) 


19.978 


20.803 


14.002 


16.528 


16.532 


57.5 


1524.9 


1-7 P(8) 


14.565 


14.593 


7.929 


9.677 


9.255 


23.5 


1547.6 


1-8 


R(3) 


10.831 


11.200 


7.458 


8.565 


8.528 


35.3 


1557.2 


1-8 R(6) 


10.993 


9.828 


4.920 


6.230 


5.864 


17.0 


1562.7 


1-8 


P(5) 


13.631 


13.810 


9.112 


10.847 


10.845 


37.2 


1581.0 


1-8 P(8) 


10.123 


9.813 


5.191 


6.384 


6.088 


17.5 


1603.5 


1-9 


R(3) 


3.851 


3.685 


2.342 


2.774 


2.773 


11.2 


1611.3 


1-9 R(6) 


3.608 


3.098 


1.520 


1.941 


1.816 




1618.2 


1-9 


P(5) 


4.742 


4.487 


2.843 


3.482 


3.459 


11.6 


1634.2 


1-9 P(8) 


3.653 


3.256 


1.630 


2.048 


1.935 


5.5 


1656.2 


1-10 


R(3) 


0.443 


0.402 


0.250 


0.312 


0.307 




1661.5 


1-10 R(6) 


0.395 


0.332 


0.161 


0.207 


0.193 




1670.0 


1-10 P(5) 


0.537 


0.487 


0.303 


0.382 


0.374 




1682.8 


1-10 P(8) 


0.425 


0.359 


0.175 


0.224 


0.209 




Pumped by 3-1 P(14) 


Pumped by 4-3 P(5) 



Table 2 — Continued 







fline 


(xl0~ 15 org/s/cm 2 ) 








Fun, 


(xl0~ 


15 erg/s/cm 2 ) 




A 


Transition 


Model A (a m 


ax) 


Model B 


Obs. a 


A 


Transition 


Model A (a„ 


,ax) 




Model B 


Obs. a 


(A) 




10/im 1mm 


10cm 


Quics. Turb. 




(A) 




10/im 1mm 


10cm 




Quics. Turb. 





1180.3 


3-1 R(12) 


2.681 


2.242 


0.848 


0.827 


0.777 




1151 


3 


4-2 R(3) 


0.262 


0.196 


0.096 


0.153 


0.170 




1214.2 


3-1 P(14) 


3.749 


2.662 


0.854 


0.919 


0.846 




1163 


8 


4-2 P(5) 


0.333 


0.244 


0.117 


0.202 


0.217 




1231.3 


3-2 R(12) 


1.839 


1.133 


0.338 


0.405 


0.354 




1202 





4-3 R(3) 


2.094 


1.680 


0.823 


1.218 


1.354 




1266.1 


3-2 P(14) 


1.651 


1.064 


0.324 


0.367 


0.331 




1214 


9 


4-3 P(5) 


2.699 


2.108 


1.025 


1.666 


1.815 




1283.1 


3-3 R(12) 


0.345 


0.203 


0.059 


0.075 


0.063 




1253 


8 


4-4 R(3) 


3.058 


2.094 


0.945 


1.597 


1.745 


6.0 


1318.7 


3-3 P(14) 


0.328 


0.195 


0.057 


0.071 


0.061 




1267 





4-4 P(5) 


3.771 


2.551 


1.161 


2.238 


2.324 


12.9 


1335.4 


3-4 R(12) 


2.896 


1.914 


0.590 


0.671 


0.605 




1306 


4 


4-5 R(3) 


0.142 


0.091 


0.041 


0.086 


0.085 




1371.5 


3-4 P(14) 


2.432 


1.746 


0.564 


0.605 


0.557 


4.0 


1319 


8 


4-5 P(5) 


0.171 


0.109 


0.049 


0.105 


0.102 




1387.7 


3-5 R(12) 


1.246 


0.772 


0.231 


0.275 


0.241 


3.4 


1359 


3 


4-6 R(3) 


3.399 


2.430 


1.105 


1.751 


1.938 


6.5 


1423.9 


3-5 P(14) 


1.121 


0.729 


0.223 


0.251 


0.227 


1.6 


1372 


9 


4-6 P(5) 


4.266 


2.993 


1.380 


2.534 


2.689 


7.2 


1439.4 


3-6 R(12) 


0.597 


0.360 


0.106 


0.130 


0.112 




1412 





4-7 R(3) 


0.802 


0.526 


0.236 


0.443 


0.466 




1475.2 


3-6 P(14) 


0.557 


0.346 


0.103 


0.121 


0.107 




1425 


6 


4-7 P(5) 


0.976 


0.637 


0.287 


0.587 


0.590 




1489.4 


3-7 R(12) 


2.669 


1.871 


0.595 


0.657 


0.600 


6.7 


1463 


8 


4-8 R(3) 


1.957 


1.360 


0.615 


1.014 


1.115 


2.3 


1524.3 


3-7 P(14) 


2.166 


1.680 


0.573 


0.599 


0.555 


5.1 


1477 


3 


4-8 P(5) 


2.440 


1.671 


0.764 


1.446 


1.515 


5.8 


1536.8 


3-8 R(12) 


0.658 


0.403 


0.120 


0.144 


0.126 




1513 


7 


4-9 R(3) 


2.463 


1.787 


0.815 


1.269 


1.408 


4.3 


1569.8 


3-8 P(14) 


0.612 


0.390 


0.118 


0.135 


0.121 




1526 


8 


4-9 P(5) 


3.125 


2.221 


1.028 


1.861 


1.986 


6.1 


1580.0 


3-9 R(12) 


1.616 


1.082 


0.335 


0.379 


0.343 


8.9 


1560 


6 


4-10 R(3) 


0.819 


0.546 


0.245 


0.440 


0.472 




1609.9 


3-9 P(14) 


1.385 


1.012 


0.331 


0.352 


0.325 


8.5 


1572 


9 


4-10 P(5) 


1.011 


0.668 


0.302 


0.604 


0.615 


6.8 


1617.0 


3-10 R(12) 


3.614 


3.050 


1.108 


1.154 


1.075 


6.8 


1602 


8 


4-11 R(3) 


4.721 


4.325 


2.165 


2.803 


3.057 


17.6 


1642.1 


3-10 P(14) 


2.722 


2.538 


1.069 


1.068 


1.001 




1613 


9 


4-11 P(5) 


6.466 


5.715 


2.925 


4.464 


4.952 


20.4 



'FUSE and HST observations by Herczeg et al. (2002). 



